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We give an overview of the present status of knowledge of the production of J/f/', '/'' and 
T in high-energy hadron collisions. We first present two early models, namely the Colour- 
Singlet Model (CSM) and the Colour-Evaporation Model (CEM). The first is the natural 
application of pQCD to quarkonium production and has been shown to fail dramatically 
to describe experimental data, the second is its phenomenological counterpart and was 
introduced in the spirit of the quark-hadron duality in the late seventies. Then, we 
expose the most recent experimental measurements of J/'i/'i ''!>' and T prompt and direct 
production at nonzero PT from two high-energy hadron colliders, the Tevatron and RHIC. 
In a third part, we review six contemporary models describing J/"!/), V' and T production 
at nonzero pT- 
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1. History: from a revolution to an anomaly 
1.1. J/ij} and the November revolution 

The era of quarkonia has started at the simultaneous discovery'^ of the J /ip 
in November 1974 by Ting et aZP at BNL and by Richter et aM at SLAC. 
Richter's experiment used the electron-positron storage ring SPEAR, whose center- 
of-momentum energy could be tuned at the desired value. With the Mark I detector, 
they discovered a sharp enhancement of the production cross section in different 
channels: e+e", /i^/i^, tt+tt^,. . . On the other hand. Ting's experiment was based 
on the high-intensity proton beams of the Alternating Gradient Synchrotron (AGS) 
working at the energy of 30 GeV, which bombarded a fixed target with the conse- 
quence of producing showers of particles detectable by the appropriate apparatus. 

In the following weeks, the Frascati group (Bacci et azP|) confirmed the presence 
of this new particle whose mass was approximately 3.1 GeV. The confirmation was 
so fast that it was actually published in the same issue of Physical Review Letters, 
ie. vol. 33, no. 23, issued the second of December 1974. In the meantime, Richter's 
group discovered another resonant state with a slightly higher mass, which was 
called*^ i)'. 

It was also promptly established that the quantum numbers of the J /^p were the 
same as those of the photon , i.e. 1 . Moreover, since the following ratio 

cross section for e+e^ — > hadrons , ^ 

R = ; I _ T — ^~ (1) 

cross section for e^e —t ji^n 

was much larger on-resonance than off, it was then clear that the J/ip did have 
direct hadronic decays. The same conclusion held for the ijj' as well. The study 
of multiplicity in pion decays indicated that ip decays were restricted by G-parity 
conservation, holding only for hadrons. Consequently, J/ip and ip' were rapidly 
considered as hadrons of isospin and G-parity -1. 

Particles with charge conjugation G different from -1 were found later. Indeed, 
they were only produced by decay of ip" and the detection of the radiated photon 
during the (electromagnetic) decay was then required. The first to achieve this task 
and discover a new state was the DASP collaboratiorP based at DESY (Deutsches 
Elektronen-Synchrotron), Hamburg, working at an e+e" storage ring called DORIS. 
This new particle, named'^ Pc, had a mass of approximately 3.5 GeV. At SPEAR, 
other resonances at 3.415, 3.45 and 3.55 GeV were discovered, the 3.5 GeV state 
was confirmed. Later, these states were shown to be G = -|-1. 

Coming back to the ratio R, it is instructive to analyse the implication in the 
framework of the quark model postulated by Gell-Mann and Zwcig in 1963. In 

''Nobel prize of 1976. 

'^We shall also make us of the name ip{2S). 

•^The name chosen here reveals that physicists already had at this time an interpretation of this 
state as bound-state of quarks c. Indeed, the symbol P follows from the classical spectroscopic 
notation for a. I = 1 state, where I stands for the angular momentum quantum number. 
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Fig. 1. Experimental status of R as of July 1974 (from Ref. 151 . 

1974 at the London conference, Richter presentea^ the experimental situation as 
in Fig.HI 

In the framework of the GeU-Mann-Zweig quark model with three quarks, a 
plateau was expected with a value of 2/3 or 2 if the quark were considered as 
"coloured" - a new concept introduced recently then The sign of a new quark - as 
the charm quark first proposed by Bjorken and Glashow in 196^ would have been 
another plateau from a certain energy (roughly its mass) with a height depending 
on its charge. Retrospectively, one cannot blame anyone for not having seen this 
"plateau" on this kind of plots. Quoting Richter, "the situation that prevailed in 
the Summer of 1974" was "vast confusion" . 

This charm quark was also required by the mechanism of Glashow, Iliopoulos and 
Maiani (GIM)2i, in order to cancel the anomaly in weak decays. The charm quark 
was then expected to exist and to have an electric charge 2/3. R was therefore to be 
10/3 in the coloured quark model, still not obvious in Fig.^ This explains why the 
discovery of such sharp and easily interpreted resonances in November 1974 was a 
revolution in particle physics. 

It became quite rapidly obvious that the J/ijj was the lowest-mass cc system with 
the same quantum numbers as photons - explaining why it was so much produced 
compared to some other members of its family. These cc bound states were named 
"charmonium" , firstly by Appelquist, De Riijula, Politzer and Glasho'vJ^in analogy 
with positronium, whose boimd-state level structure was similar. 

At that time, the charm had nevertheless always been found hidden, that is in 
charm-anti-charm bound states. In order to study these explicitly charmed mesons, 
named D, the investigations were based on the assumption that the D was to 
decay weakly and preferentially into strange quarks. The weak character of the 
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decay motivated physicists to search for parity-violation signals. The first resonance 
attributed to meson was found in if+Tr"*" decay by Goldhaber et aZ0in 1976. A 
little later, _D+ and were also discovered as well as an excited state, Z?*, with 
a mass compatible with the decay ip'" ~^ D^D* . And, finally, the most conclusive 
evidence for the discovery of charmed meson was the observation of parity violation 
m D decayOT 

To complete the picture of the charm family, the first charmed 
baryon was discovered during the same yeaJl^l The quarks were not anymore just 
a way to interpret symmetry in masses and spins of particles, they had acquired an 
existence. 



1.2. The bottomonium family 

In the meantime, in 1975, another brand new particle was unmasked at SLAG by 
Perl et a^^. This was the first particle of a third generation of quarks and leptons, 
the r, a lepton. Following the standard model, more and more trusted, two other 
quarks were expected to exist. Their discovery would then be the very proof of 
the theory that was developed since the sixties. Two names for the fifth quark were 
already chosen: "beauty" and "bottom" , and in both cases represented by the letter 
h\ the sixth quark was as well already christened with the letter t for the suggested 
names "true" or "top". 

The wait was not long: two years. After a false discovery'^ of a resonance at 6.0 
GeV, a new dimuon resonance similar to J/V' and called T was brought to ligh t at 
Fermilab, thanks to the FNAL proton synchrotron accelerator, by Herb et a/El 
with a mass of 9.0 GeV; as for charmonia, the first radial excited state (T(25)) 
was directly founcfl^ thereafter. Again, the discovery of a new quarkonium was a 
decisive step forward towards the comprehension of particle physics. 

Various co nfirm ations of these discoveries were not long to come. The 35 state 
was then founcJl^ at Fermilab as well as an evidence that the AS state was lying 
above the threshold for the production of B mesons. The latter was confirmed at the 
Gornell e+e~ storage ring with the CLEG detector. The first evidence for B meson 
and, thus, for unhidden h quark, was also brought by the GLEO coUaboratiorP-^ in 
1980. One year later, "the first evidence for baryons with naked beauty" (sic) was 
reported by GERN physicist^^. 

Another decade was needed for the discovery of the sixth quark which was 
definitely christened "top" . Indeed, in 1994, the GDF GoUaboration found the first 
evidence for it at the Tev atron collider at Fermilabi^. The discovery was published 
m 1995 both by GDjEHI 

and Unfortunately, due to its very short hfetime, 

this quark cannot bind with its antiquark to form the toponium. To conclude this 
historical prelude, we give the spectra (Figs. 2 & 3) of the cc and bb systems as well 
as two tables (Tables 1 & 2) summing up the characteristics of the observed states 

"^Nonetheless, this papei^^l first suggested the notation T for any "onset of high-mass dilepton 
physics" . 
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as of today. 



Table 1. Properties of charmonia (cf. Ref. 1221 . 



Meson 






Mass (GeV) 


Tmm (keV) 


Vc 


1 ^So 


0-+ 


2.980 


N/A 




1 ^Si 


1 — 


3.097 


5.40 


XcO, Xcl, Xc2 


1 ''Po,i,2 


0++,l++,2++ 


3.415,3.511,3.556 


N/A 


he 


1 


1+- 


3.523 


N/A 


Vc{2S) 


2 


0-+ 


3.594 


N/A 


■ip' 


2 35i 


1 — 


3.686 


2.12 



V(2S) 




hadrons hadrons y* radiative 



jPC = 0-+ 1— 0++ 1++ 1+- 2++ 

Fig. 2. Spectrum and transitions of the charmonium family (Reprinted figure from Ref. 1221 with 
permission of Elsevier. Copyright (2004).). 



Table 2. Properties of bottomonia (cf. Ref. 221 . 



Meson 






Mass (GeV) 


Tmm (keV) 


T(15) 


1^5i 


1 — 


9.460 


1.26 


X60: X61 > X62(l-P) 


1 ^ft,l,2 


0++,l++,2++ 


9.860,9.893,9.913 


N/A 


T(25) 


2 


1 — 


10.023 


0.32 


X60. Xbi, Xb2 (2P) 




0++,l++,2++ 


10.232,10.255,10.269 


N/A 


T(35) 


3 ^Si 


1 — 


10.355 


0.48 



1.3. Early predictions for quarkonium production 
1.3.1. The Colour-Singlet Model 

This model is° the most natural appHcation of QCD to heavy-quarkonium produc- 
tion in the high-energy regime. It takes its inspiration in the factorisation theorem of 



•^before the inclusion of fragmentation contributions. 
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Fig. 3. Spectrum and transitions of the bottomonium family (Reprinted figure from Ref.ra with 
permission of Elsevier. Copyright (2004).). 



where the hard part is calculated by the strict application pQCD 
and the soft part is factorised in a universal wave function. This model is meant to 
describe the production not only of J/'0, ip{2S), T(IS'), T(25') and T(3S'), i.e. the 
'^Si states, but also the singlet S states rjc and rjh as well as the P (x) and D states. 

Its greatest quality resides in its predictive power as the only input, apart from 
the PDF, namely the wave function, can be determined from data on decay pro- 
cesses or by application of potentia l models. Nothing more is required. It was first 
applied to hadron colliderj2321ESI^ then to electron-proton coUiderJSSI. The cross 
sections for '^S'l states were then calculated, as well as also for 77 and Xi for charmo- 
nium and for bottomonium. These calculations were compared to ISR and FNAL 
data from = 27 GeV to y/s = 63 GeV fo r which the data extended to 6 GeV 
for the transverse momentum. UpdateJ^^ of the model to describe collisions at 
the CERN pp collider {y/s — 630 GeV) were then presented. At that energy, the 
possibility that the charmonium be produced from the decay of a beauty had ron 
was becoming competitive. Predictions for Tevatron energies were also madJ^. 

In order to introduce the reader to several concepts and quantities that will be 
useful throughout this Review, let us proceed with a detailed description of this 
model. 



proven for some definite cases, e.g. Drell-Yan process. 
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1.3.2. The Model as of early 90's 

It was then based on several approximations or postulates: 

• If we decompose the quarkonium production in two steps, first the creation 
of two on-shell heavy quarks [Q and Q) and then their binding to make the 
meson, one postulates the factorisation of these two processes. 

• As the scale of the first process is approximately + Pr, one considers 
it as a perturbative one. One supposes that its cross section be computable 
with Feynman-diagram methods. 

• As we consider only bound states of heavy quarks (charm and bottom 
quarks), their velocity in the meson must be small. One therefore supposes 
that the meson be created with its 2 constituent quarks at rest in the meson 
frame. This is the static approximation. 

• One finally assumes that the colour and the spin of the QQ pair do not 
change during the binding. Besides, as physical states are colourless, one 
requires the pair be produced in a colour- singlet state. This explains the 
name Colour-Singlet Model (GSM). 

In high-energy hadronic collisions, the leading contribution comes from a gluon 
fusion process; as the energy of the collider increases, the initial parton momentum 
fraction Xi needed to produce the quarkonium decreases to reach the region in x 
where the number of gluons becomes much larger than the number of quarks. One 
has then only six Feynman diagrams for the '^Si states production associated with 
a gluon^ (see Fig.0J. 



3i 



Si . 





Fig. 4. The 6 diagrams for gg ^Sig at LO within the CSM. 



One usually starts with A4 (p) , the perturbative amplitude to produce the heavy- 
quark pair on-shell with relative momentum p and in a configuration similar to 
the one of the meson. To realise the latter constraint, one introduces a projection 

^This is the dominant process when the transverse momentum of the meson is non- vanishing. 
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operator'^; the amplitude M{p) is then simply calculated with the usual Feynman 
rules. 

The amplitude to produce the meson is thence given by 



where is the usual Schrodinger wave-function. 

Fortunately, one does not have to carry out the integration thanks to the the 
static approximation which amounts to considering the first non-vanishing term of 
A when the perturbative part M. is expanded in p. For S'-wave, this gives 



where ^ is the wave-function in coordinate space, and ^'I^^q or ^'(0) is its value at 
the origin. For P-waves, ^(0) is zero, and the second term in the Taylor expansion 
must be considered; this makes appear \l/'(0). ^'(0) (or ^''(0)) is the non-perturbative 
input, which is also present in the leptonic decay width from which it can be ex- 
tracted. 

If the perturbative part M (p) is calculated at th e leading or der in as , this will 
be referred to as the Leading Order GSM (LO CSMj^ZEHEHEDl. 

1.3.3. The Colour Evaporation Model 

This model was initially introduced in IQIll ^^ l ^ ^ Q^^j^ ^g^g revived in 1996 by 
Halzen et aZ.^"^ -SIcontrarily to the CSM, the heavy-quark pair produced by the 
perturbative interaction is not assumed to be in a colour-singlet state. One simply 
considers that the colour and the spin of the asymptotic QQ state is randomised 
by numerous soft interactions occurring after its production, and that, as a conse- 
quence, it is not correlated with the quantum numbers of the pair right after its 
production. 

A first outcome of this statement is that the production of a ^5*1 state by one 
gluon is possible, whereas in the CSM it was forbidden solely by colour conservation. 
In addition, the probability that the QQ pair eventually be in a colour-singlet state 
is therefore ^, which gives the total cross section to produce a quarkonium: 



This amounts to integrating the cross section of production of QQ from the threshold 
2toq up to the threshold to produce two charm or beauty mesons (designated here 
by qQ) and one divides by 9 to get the probability of having a colour-singlet state. 

''In fact, this amounts to associate a 7® matrix to pseudoscalars, 7*^ to vectors, etc. 




(2) 
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The procedure to get the cross section for a definite state, for instance a J/'0, 
is tantamount to "distributing" the cross sections among aU states: 

(5) 

The natural value for pj/^, as well as for the other states in that approximation, is 
the inverse of the number of quarkonium lying between 2toc and 2m jj. This can be 
refined by factors arising from spin multiplicity arguments or by the consideration 
of the mass difference between the produced and the final states. These are included 
in the Soft-Colour-Interactions approach (SCI) (see section |OJ- 

By construction, this model is unable to give information about the polari- 
sation of the quarkonium produced, which is a key test for the other mo del 
Furthermore, nothing but fits can determine the values to input for p. Considering 
production ratios for charmonium states within the simplest approach for spin, we 
should have for instance cr[ric] : a[J/Tp] = 1:3 and cr[xco] ■ f [Xci] : <^[Xc2] = 1:3:5, 
whereas deviations from the predicted ratio for Xci and Xc2 have been observed. 
Moreover, it is unable to describe the observed variation - from one process to 
another - of the production ratios for charmonium states. For example, the ratio 
of the cross sections for Xc and J/iJj differs significantly in photoproduction and 
hadroproduction, whereas for the CEM these number are strictly constant. 

All these arguments make us think that despite its simplicity and its phenomeno- 
logical reasonable grounds, this model is less reliable than the CSM. It is also instruc- 
tive to point out that the invocation of reinteractions after the QQ pair production 
contradicts factorisation, which is albeit required when the PDF are used. However, 
as we shall see now, the CSM has undergone in the nineties a lashing denial from 
the data. 



1.4. xp' anomaly 

In 1984, Halzen et aZ.I^ noticed that charmonium production from a weak B decay 
could be significant at high energy - their prediction were made for ^/s = 540 GeV 
- and could even dominate at high enough p^. This can be easily understood having 
in mind that the B meson is produced by fragmentation of a & quark - the latter 
dresses with a light quark -. And to produce only one b quark at high px is not so 
burdensome; the price to pay is only to put one quark propagator off-shell, instead 
of two for gg —^ ipg. 

This idea was confirmed by the calculations of Glover et al\^, which were used 
as a benchmark by the UAl CoUaboratiorP^. After the introduction of a if factor 
of 2, the measurements agreed with the predictions; however the px slope was not 
compatible with 6-quark fragmentation simulations. 

From 1992, the CDF collaboration undertook an analysis of?/)' production. They 
managed to extract unambiguously the prompt component of the signal (not from 



'In the following, ip stands for both J/4> and ip' . 



10 Jean- Philippe Lansberg 



B decay) using a Silicon Vertex Detector 

The preliminary'^ results*^^ showed an unexpectedly large prompt component. 
For the ijj' , the prompt cross section was orders of magnitude above the predictions 
of the LO CSM (compare the data to the dashed curve on Fig. (right)) . This 
problem was then referred to as the tp' anomaly. For J/ip, the discrepancy was 
smaller (Fig. (left)), but it was conceivably blurred by the fact that a significant 
part of the production was believed to come from Xc radiative feed-down, but no 
experimental results were there to confirm this hypothesis. 




Fig. 5. Comparison between preliminary measurements of Ref. 1401 from CDF and the cross 
sections obtained by Braaten et al. for LO CSM (dashed curves) of J/'^p (left) and -)/)' (right), as 
well as for CSM fragmentation contribution (solid curves). In each case, the two curves depict the 
extremum values obtained by varying parameters such as rric and the different scales: fiji, fip, 
t^frag (Reprinted figures from Ref. I42l with permission of Elsevier. Copyright (1994).) 



1.4.1. Fragmentation within the CSM: the anomaly confirmed 

The prediction from the LO CSM for prompt ip being sig nificantly below these pre- 
liminary measurements by CDF, Braaten and Yuan pointed out in 1993 that 
giuon fragmentation processes, even though of higher order in as, were to prevail 
over the LO CSM for S'-wave mesons at large pr, in same spirit as the produc- 
tion from B was dominant as it came from a fragmentation process. Following this 
paper, with Cheung and Flemin^Sl they considered the fragmentation of c quark 
into & ij) in Zq decay. From this calcu lati on, they extracted the corresponding frag- 
mentation function. In another papen^l they considered gluon fragmentation into 
P-wave mesons. All the tools were then at hand for a full prediction of the prompt 
component of the J I iJj an d "0' at the Tevatron. This was realised simultaneously 
by Cacciari and Greco EHl and by Braaten et a/PS Let us now review briefly the 
approach followed. 

■■The details of the analysis are given in the following section. 
The final results -confirming the preliminary ones- were in fact published in 1997 El 
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To all orders in ofg, we have the following fragmentation cross section for a 
quarkonium Q: 

^ /■! p 

'^Q{P)-y2 dzdai{ — ,^frag)Di^Q{z,flfrag)- (6) 

The fragmentation scale, fJ-frag, is as usual chosen to avoid large logarithms of 
Pr/fJ-frag in <^i{^7 fJ-frag): that is fJ-frag — Pt ■ The Summation of the corresponding 
large logarithms of ^j. frag /mo appearing in the fragmentation function is realised 
via an evolution equation^^"^. 




The interesting point raised by Braaten and Yuan is that the fragmentation 
functions can be calculated perturbatively in at the scale fJ^frag = 2toq. For 
instance, in the case of gluon fragmentation into an rjc, the trick is to note that 

dzDz(z,mc) is the ratio to the rates for the well-known gg —> gg process and 
99 ~^ Vc99 (Fig- Ell- After some manipulations, the fragmentation function can be 
obtained from this ratio by identifying the integrand in the z integral. This gives : 

Dg^,,^iz,2m,) = ia2(2m,)M^[3^_2z2 + 2(l-z)ln(l-z)]. (7) 
6 m-^ 

The other fragmentation functions of a given parton i into a given quarkonium Q, 
Di^Q{z, ^ifrag), were obtained in the same spirit. For the Tevatron, the differential 
cross section versus px of various CSM fragmentation processes are plotted in Fig. [3 
(left). 

The prompt component of the J/iJj and the direct component of the -ip' could in 
turn be obtained and compared with the preliminary data of CDF (see the solid 
curves in plots in Fig.^above). For the J/ip, the previous disagreement was reduced 
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Fig. 7. (left) Differential cross section versus of various CSM fragmentation processes for 
J/i/i to be compared with the LO contributions (Reprinted figure from Ref. H2| with permission of 
Elsevier. Copyright (1994).). (right) Differential cross section versus px of the CSM (fragmentation 
and LO) production to be compared with the direct production of J/i/i from CDF (Reprinted figure 
from Ref. ICTl with permission of Elsevier. Copyright (2001).) 



and could be accounted for by the theoretical and experimental uncertainties; on the 
other hand, for the ip', the disagreement continued to be dramatic. The situation 
would be clarified by the extraction of the direct component for J/V', for which 
theoretical uncertainties are reduced and are similar to those for the tp' . 

The CDF collaboration undertook the disentanglement of the direct J/tp 
signaPSI They searched for J/ip associated with the photon emitted during this 
radiative decay: the result was a direct cross section 30 times above expectations 
from LO CSM plus fragmentation. This was the confirmation that the CSM was 
not the suitable model for heavy-quarkonium production in high-energy hadron 
collisions. 

It is a common misconception of the CSM to believe that the well-known factor 
30 of discrepancy between data and theory for direct production of J/ip arises when 
the data ar e comp ared with the predictions for the LO CSM following Baier, Riickl, 
Chang, . . . ) ' l ^v l ^y l tuned to the right energy. As you can see on Fig. |3 (right), the 
factor would be rather of two orders of magnitude at large pt for J/iIj. The same 
conclusion holds also for ip' (see Fig. (right)). 

It is also worth pointing that, in the CSM, the direct component of tp {^Si 
charmonia) produced by fragmentation is mainly from c-quark fragmentation (see 
Fig. 13 (left)) as soon as Pt reaches 5 GeV and the ag penalty of the gluon frag- 
mentation is not compensated anymore by the c-quark mass. It was further pointed 
out in 2003 by Qiao''"Jthat sea-quark initiated contributions could dominate in the 
fragmentation region (large p-r)- 
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2. Review of contemporary measurements of direct production of 
J/'0, xp' and T from the Tevatron and RHIC 

2.1. Foreword 

Limiting ourselves to high-energy coUisions, the most recent published results for 
quarkonium hadroproduction come from two accelerators: 

(1) The Tevatron at Fermilab which - as stated by its name - runs at TeV energy 
with proton-antiproton collisions. For Run I ("Run lA" 1993-94 and "Run IB" 
1995-96), the energy in cm. was 1.8 TeV. For Run II, it has been increased to 
1.96 TeV. The experimental analyses for this Run arc still being carried out. 

(2) RHIC at BNL running at 200 GeV for the J/V' study with proton-proton colli- 
sions. 

It is worth pointing out here that high-energy p-p and p-p collisions give similar 
results for the same kinematics, due to the small contribution of valence quarks. 



2.2. Different types of production: prompt, non-prompt and direct 

As we have already explained, the detection of quarkonia proceeds via the iden- 
tification of their Icptonic-decay products. We give in Table. 01 the relative decay 
widths into muons. 



Table 3. Table of branching 
ratios in dimuons (Ref. 1221 . 



Meson 


T{^l+ ^l-)/T{total) 


J/t/. 


0.0588 ±0.0010 


V.(25) 


0.0073 ± 0.0008 


T(1S) 


0.0248 ± 0.0006 


T(2S) 


0.0131 ± 0.0021 


T(3S) 


0.0181 ± 0.0017 



Briefly, the problem of direct J /ip production separation comprises three steps: 

• muon detection; 

• elimination of J /tp produced by hadrons containing b quarks ; 

• elimination of Xc radiative-decay production; 

Table 4 summarises the different processes to be discussed and the quantities 
linked. 
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Table 4. Different processes involved in J/i/i production accompanied by quantities used in 
the following discussion. 



1^* step 


2"'^ step 


3'''' step 


Type 


Associated 
quantity 


pp — > cc + Jf 


cc — » J/4> 
cc Xc 

cc — > ip' 


Xc J/tp + 7 


Direct prod. 
Prompt prod, by 

decay ofxc 
Prompt prod, by 

decay of i/>' 




pp — > 6c + X 
be ^ ce + £~ + 

etc . . . 


cc ^ J/l/l 
cc ^ Xc 


Xc J/V' + 7 




MorFt) 



Let us explain the different fractions that appear in the table: 

• F{p,^,^p{2S))"'^^ is the prompt fraction of J/tp that do not come x, neither 
from -0(25*), i.e. the direct fraction. 

• F[fiYJ^g.^ is the prompt fraction of J/ij} that come from ip(2S). 

• F{fi)'lJ^ is the prompt fraction of J /ijj that come from x- 

• Ff, (or /f,) is the non-prompt fraction or equally the fraction that come for h 
quarks. 

Concerning tp{2S), due to the absence of stable higher excited states likely to 
decay into it, we have the summary shown in Table 5 for the different processes to 
be discussed and the quantities linked: 



Table 5. Different processes involved in i/i' production accompanied by quantities used in 
the following discussion. 



1"* step 


2"'' step 


3'''* step 


Type 


Associated 
quantity 


pp ^ cc + X 


cc ^ V(25) 




Direct/Prompt prod. 




pp ^ be + X 
be CC + + vi 
etc . . . 


cc V(25) 






MorFt) 



Let us explain the different fractions that appear in the corresponding table: 

• F{fi)^'''^^^ is the prompt fraction of ip{2S), i.e. the direct production. 

• Fb (or fh) is the non-prompt fraction or equally the fraction that come for b 
quarks. 

As can be seen in Table. El the leptonic branching ratio of T are also relatively 
high. The detection and the analysis of the bottomonia are therefore carried out 
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in the same fashion. For the extraction of the direct production, the 6-quark feed- 
down is obviously not relevant, only the decays from stable higher resonances of the 
family are to be considered. 

All the quantities useful for the bottomonium discussion are summarised in Table 

6: 



Table 6. Different processes involved in T{n5) in = 1, 2, 3) production accompanied by quan- 
tities used in the following discussion. 



1"* step 


2"'* step 


step 


Type 


Associated 
quantity 


pp^hh + X 


bb r{nS) 
bh Xb 

bb r{n's) 


Xb T{n5) + 7 
T{n5) -> r{n'S) + X 


Direct prod. 
Prod, by 

decay of Xb 
Prod, by 
decay of 

r{n's) 


~t(nS) 
direct 
„T(nS) 

„T(nS) 



Let us explain the different fractions that appear in the latter table: 

• is the direct fraction of T{nS). 

• ^Xi,^"^^ is the fraction of T{nS) that come from Xb- 

• i^x/n'S) i^ fraction of T(jiS) that come from a higher T(n'S). 

2.3. CDF analysis for ip production cross sections 

The sample of collisions amounts to 17.8±0.6 pb^^ at yG ~ 1.8 TeV^5 For J/ip, 
the considered sample consists however of 15.4±0.6 pb^^ of integrated luminosity^. 

The number of -0 candidates is therefore determined by fitting the mass dis- 
tribution of the muons in the cm. frame after subtraction of the noise. The mass 
distribution is fit to the signal shape fixed by simulation and to a linear background. 
The fit also yields the mass of the particle and a background estimate. 

From the branching ratio and after corrections due to the experimental effi- 
ciency, the number of produced particles is easily obtained from the number of the 
candidates. 

For the present study of CDF, the fits are reasonably good for each p^-bin, the 
per degree of freedom ranging from 0.5 to 1.5. The measured width of the mass 
peak was from 17 MeV to 35 MeV for pr from 5 to 20 GeV. 

Approximatively, 22100 J/tp candidates and 800 ip' candidates above a back- 
ground of 1000 events are observed. The J/ip efficiency is 97.0 ± 0.2% and the ip' 
one is 92.3 ± 0.2#2. 

''This difference is due to the subtraction of data taken during a period of reduced level 3 tracking 
efficiency. These data have however been taken into account for '4>{2S) after a correction derived 
from the J/4> sample. 



16 Jean- Philippe Lansberg 



The integrated cross sections are measured to be 

ct( J/V') • B{J/ip fi+fi-) = 17.4 ± 0.l{stat.)+l-^g{syst.) nb (8) 

(7{^{2S)) ■ B{ij(2S) ^ n+ir) = 0.57 ± QM{stat.)tl-ll{syst.) nb. (9) 

where cr(V') = (^{pp — > iIjX,pt{iP) > 5 GeV, |77(V')| < 0.6). 

2.3.1. Disentangling prompt charmonia 

As already seen, prompt ip^s are the ones which do not come from the decay of 
B mesons. Their production pattern has the distinctive feature, compared to non- 
prompt ones, that there exists a measurable distance between the B production 
vertex and its decay into charmonium. 

To proceed, CDF uses the SVX, whose resolution is 40/im, whereas B lifetime is 
ctb ~ 450/^m. Muons are constrained to come from the same point which is called 
the secondary vertex, as opposed to the primary vertex, that is the collision point 
of protons. Then the projection of the distance between these two vertices on the 
ip momentum, L^y, can be evaluated. It is converted into a proper time equivalent 
quantity using the formula, cr — p^^^""" , where Fcorr is a correction factor, 

estimated by Monte-Carlo simulations, which links the ^ boost factor PtI to that 

of the sEa 

The prompt component of the signal is parametrised by ct = (a single ver- 
tex), the component coming from B decay is represented by an exponential, whose 
lifetime is Tb and which is convoluted with the resolution function. 

The CT distribution is fit in each p^-bin with an unbinned log-likelihood function. 
The noise is allowed to vary within the normalisation uncertainty extracted from the 
sidebands. The fraction of coming from b, fbipr), obtained by CDF is displayed 
as a function of pr in Fig. |H| 

The ip production cross section from B decays is thus extracted by multiplying*^ 
fl^^ipr) by the inclusive tp cross section. Multiplying the latter by (1 — //**(pt)), 
one obviously gets the prompt-production cross section (cf. Fig. EJ. 

We remind the reader than for ip{2S) the prompt production identifies with the 
direct one. 

2.3.2. Disentangling the direct production ofj/ip 

The problem here is to subtract the J/tp coming from Xc decay, assuming that this 
is the only source of prompt J/^ besides the direct production after subtraction of 



''In order to reduce statistical fluctuations, is fit by a parabola weighted by the observed shape 
of the cross sectiorPSI. 



J/i/i, V' '"^'i ^ Production at Hadron Colliders: a review 17 




° 6 8 10 12 14 16 18 20 
Pt(iA) [GeV/d 

Fig. 8. Fraction of i/i from iJ decay as a function of (Reprinted figure from Ref. 1411 with 
permission of American Physical Society. Copyright (1997)). 
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B from the prompt component of -0 (data from Ref. 1411 . 



F{p)J^2S)' prompt fraction of J/tp that come from tp{2S). The latter is evaluated 
by CDF from the ip{2S) cross section from the previous section and from Monte- 
Carlo simulation of the decays ip{2S) — > J/tpX where X — tttt, t] and 7r°. The 
delicate point here is the detection of the photon emitted during the radiative decay 
of the Xc- 

The sample they use is 34367 J/ip from which 32642 ± 185 is the number of real 
Jf-ip when the estimated background is removed. 

The requirements to select the photon were as follow: 
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• an energy deposition of at least 1 GeV in a cell of the central electromagnetic 
calorimeter; 

• a signal in the fiducial volume of the strip chambers (CES); 

• the absence of charged particles pointing to the photon-candidate cell (the no- 
track cut). 

The direction of the photon is determined from the location of the signal in the 
strip chambers and from the event interaction point. All combinations of the J/^p 
with all photon candidates that have passed these tests are made and the invariant- 
mass difference defined as AM = M{^i+fi~-/) - M{fj.+^j,-) can then be evaluated. 
As expected, the distribution AM shows of a clear peak from Xc decays is visible 
near AM — 400 MeV. Yet, distinct signals for Xci and Xc.2 are not resolved as the 
two states are separated by 45.6 MeV and as the mass resolution of the detector is 
predicted to be respectively 50 and 55 MeV. 

Eventually, the AM distribution obtained from the data is fit with a gaussian 
and with the background fixed by the procedure explained above but with a free 
normalisation. The parameter of the gaussian then leads to the number of signal 
events: 1230 ± 72 Xc- 

The analysis of the direct J/ip signal is done within four pT-hins: 4 < p':^^^ < 6, 
6 < p^/''' < 8 , 8 < p^/''' < 10 and p^^''' > 10 GeV. For p^^'^ > 4.0 GeV and 
l^^^/V"! < 0.6, CDF finds that the fraction of J/tp from Xc is then 

pJ/4: ^ 27.4% ± 1.6%istat.) ± 5.2%{syst.) (10) 

The last step now is the disentanglement of the prompt Xc production, that is the 
determination of F{p)'^^^ . Let us here draw the reader's attention to the fact that 
by selecting prompt J/ ^p (cf. section 12.3. l|l , J/-0 produced by non-prompt Xc have 
also been eliminated; it is thus necessary to remove only the prompt production by 
Xc decay, and nothing else. This necessitates the knowledge of F{p)x^ , the prompt 
fraction of J/^ that come from x- 

The latter is calculated as follows: 

Fipy^/^ = F^/^'^^, (11) 

where N^, TV/^'^ are the number of reconstructed Xc and J/ip from 6's, F^, F^^^ 
arc the corresponding fractions. 

F^^^ (or /;,) is known as seen in the section 1^.3. II F^ is obtained in same way 
and is 17.8% ± 0.45% for pr > 4.0 GeV. Consequently, 

Fip):!^/'^ = 29.7% ± 1.7%{stat.) ± b.7%{syst.) (12) 

and its evolution as a function oi pt is shown in Fig. ^| (left) . It is also found that 
the fraction of directly produced J/ip is 

^^i(tt = 64%±6%, (13) 
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and is almost constant from 5 to 18 GeV in px (see Fig. ^| (left)). We therefore 
conclude from the analysis of CDF that the direct production is the principal con- 
tribution to J /tp. 
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Fig. 10. (left) Fractions of J/i/) with the contribution of fe's removed. The error bars correspond 
to statistical uncertainty. The dashed lines show the upper and lower bounds corresponding to 
the statistical and systematic uncertainties combined (from 15^ . (right) Differential cross section 
for prompt production of J/rp ^ M^M^ ^ function of Pt (Reprinted figure from Ref . 1491 with 
permission of American Physical Society. Copyright (1997)). 



In order to get the cross section of direct J/i/j production, it is sufficient now to 
extract the contribution of ip' obtained by Monte-Carlo simulation and of Xc ob- 
tained by multiplying the cross section of prompt production by the factor F{p)'^^^ , 
which is a function of p'^^^ . The different cross sections are displayed in Fig. [TUI 
(right). 



2.3.3. Prompt J/ip production at ^/s = 1.96 TeV 

The first results of the run II for prompt J/ijj production at ^/s = 1.96 TeV have 
recently been published in Ref. i54i They correspond to an integrated luminosity of 
39.7 pb~^. The inclusive J/ip cross section was measured for Pt from to 20 GeV 
and the prompt signal was extracted from Pt = 1.25 GeV. The rapidity domain is 
stih -0.6 <y < 0.6. 

We do not give here the details of the experimental analysis which is thoroughly 
exposed in Ref. iMl The prompt-signal extraction follows the same lines as done for 
the analysis previously exposed. The prompt J/4' cross section obtained is plotted 
in Fig. ini 
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Fig. 11. Prompt J/ip cross section as measured by CDF at y/s = 1.96 TeV (statistical and 
systematical errors have been combined quadratically for this plot); the data are from Ref . 1541 



2.4. CDF measurement of the T production cross sections 

In this section, the results by CDF on T production in pp at -^/i =1.8 TeV are 
exposed. These results were exposed in two Letters fRefs. 1551 and we shall 
mainly focus on the second one, which considered data collected in 1993-95 and 
corresponding to an integrated luminosity of 77 ±3 pb^^. The number of candidates 
are 4430 ± 95 for T(15), 1114 ± 65 for T(25) and 584 ± 53 for T(3S'). 

The cross section for T(IS') is shown in Fig. [H (left), for T{2S) in Fig. [HI 
(middle) and for T(3S') in Fig. [Tl(right). 
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Fig. 12. left: Differential cross section of T(IS') fi /x+ as a function of for \y\ < 0.4. middle: 
Idem for T{2S) fj,~ fi'^ . right: Idem for T(35) — > A'^M^ (data from Ref. 1561 . 
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2.4.1. Disentangling the direct production ofTilS) 

The analysiJS3 presented here is for the most part the same as described in sec- 
tion 12.3.21 It is based on 90 pb^^ of data coUected during the 1994-1995 run. The 
measurement has been constrained to the range pt > 8.0 GeV because the energy 
of the photon emitted during the decay of Xb decreases at low pt and ends up to 
be too small for the photon to be detected properly. In the same spirit, analysis 
relative to T(2S') has not been carried out, once again because of the lower energy 
of the radiative decay. Concerning T(3S'), except for the unobserved Xh(3P) (which 
is nevertheless supposed to be below the BB threshold), no states are supposed to 
be possible parents. 

A sample of 2186 T(IS') candidates is obtained from which 1462 ± 55 is the 
estimated number of T(15') after subtraction of the background. In the sample 
considered, photons likely to come from Xb decay are selected as for the J/ijj case, 
except for the energy deposition in the central electromagnetic calorimeter, which 
is lowered to 0.7 GeV. 

For pj'^'^^ > 8.0 GeV and Iry^^^^^l < 0.7, the fractions of T(IS') from Xb{lP) 
and from Xf,(2P) are measured by CDF to be 

p^(^^\ =27.1% ± 6.9%(stat.) ± 4.4%(s?/si.), 
XbK^t-, ^^^^ 

F^^^i^P) =10.5% ± AA%{stat.) ± lA%{syst.). 

The feed-down from the S'-waves T(2S') and T(3S') is obtained by Monte-Carlo 
simulations of these decays normalised to the production cross sections discussed in 
sectionE31 It is found that for p^^-^^'' > 8.0 GeV the fraction of T(IS') from T(2S') 
and T{3S) are respectively 

^T{2S) ~ 4 8 '^J -^T(3S)^ = 0.8%j^Q'4%. (15) 

Concerning the unobserved Xbi'^P), a maximal additional contribution is taken 
into account by supposing than all the T(3S') are due to Xf,(3i-') and from theoretical 
expectation for the decay of this state, a relative rate of T(IS') from Xf,(3P) can 
obtained. This rate is less than 6%. 

Eventually the fraction of directly produced T(IS') is found to be 

PjJ^-lct = 50.9% ± 8.2%(stat.) ± 9.0%{syst.). (16) 
2.5. Polarisation study 

As the considered bound states, and T are massive spin-1 particles, they have 
three polarisations. In addition to measurements of their momentum after their 
production, the experimental set-up of CDF is sufficiently refined to provide us 
with a measurement of their spin alignment through an analysis of the angular 
distribution of the dimuon pairs from the decay. 
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The CDF collaboration has carried out two analyses, one for the V' statea^^ - 
for which the feed-down from b decay has been subtracte d, b ut not the indirect 
component in the case of the J/ip - and another for T{nS)'^. 

In the following, we shall proceed to a brief outline of the analysis and give its 
main results. 




Fig. 13. Definition of the angle 9 used in the polarisation analysis of a quarkonium Q. 



The polarisation state of the quarkonium can be deduced from the angular 
dependence of its decay into ^'^fi~. Taking the spin quantisation axis along the 
quarkonium momentum direction in the pp cm. frame, we define 9 as the angle be- 
tween the /i+ direction in the quarkonium rest frame and the quarkonium direction 
in the lab frame (see Fig. I13|l . Then the normalised angular distribution /(cos 6*) is 
givcn"^ by 

I{cos9)^—^—{l + acosH) (17) 
Z[a + 6) 

where the interesting quantity is 



(18) 

a = means that the mesons are unpolarised, a — +1 corresponds to a full trans- 
verse polarisation and a = — 1 to a longitudinal one. 

As the expected behaviour is biased by muons cuts - for instance there exists a 
severe reduction of the acceptance as 9 approaches and 180 degrees, due to the pt 
cuts on the muons -, the method followed by CDF was to compare measurements. 



•^For a derivation, see the Appendix A of R.ef . 1581 
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not with a possible (1 + a cos^ 9) distribution, but with distributions obtained 
after simulations of quarkonium decays taking account the geometric and kinematic 
acceptance of the detector as well as the reconstruction efficiency. 

The parts of the detector used are the same as before, with the additional central 
niuon upgrade (CMP) outside the CMU. 

2.5.1. Study of the tp 's polarisation by CDF 

We give here the results relative to the CDF analysis published in 2000 The 
data used correspond to an integrated luminosity of 110 pb~^ collected between 
1992 and 1995. 

Disentangling prompt production for ip The measured fraction of J/tp mesons 
which come from &-hadron decay, F^^^ , is measured to increase from (13.0 ±0.3)% 
at pT = 4 GeV to (40 ± 2)% at 20 GeV and for V'' mesons, Ff is (21 ± 2)% at 5.5 
GeV and (35 ± 4)% at 20 GeV. 

Within a 3-standard-deviation mass window around the J lijj peak, the data 
sample is of 180 000 J/V' events. In order to study the effect of pr, the data are 
divided into seven p^-bins from 4 to 20 GeV. Because the number of i}}' events is 
lower, data for are divided into three pT-bins from 5.5 to 20 GeV. 

J/t/j polarisation measurement The polarisation is obtained using a fit of 
the data to a weighted sum of transversely polarised and longitudinally polarised 
templates. The weight obtained with the fit provides us with the polarisation. Ex- 
planations relative to procedure used can be found in Refs. 1591 and 1581 

Note however that the polarisation is measured in each pT bin and that separate 
polarisation measurements for direct J/ij} production and for production via Xc and 
■0' decays was found to be unfeasible. Let us recall here that Xc and ij}' were shown 
to account for 36±6% of the prompt production (cf. Eq. l(T^ ') and to be mostly 
constant in the considered pr range. 

Except in the lowest px bins, the systematic uncertainties are much smaller 
than the statistical one. The values obtained for aprompt and cxfromB are given in 
Table. Hand plotted in Fig.m 



Table 7. Fit results for J/i/> polarisation, with statistical and systematic uncertainties 
(Ref.EHJ. 



Pt bin (GeV) 


Mean Pt (GeV) 


^Prompt 




4- 5 

5- 6 

6- 8 
8-10 
10 - 12 
12 - 15 
15 - 20 


4.5 
5.5 
6.9 
8.8 
10.8 
13.2 
16.7 


0.30 ±0.12 ±0.12 
0.01 ±0.10 ±0.07 
0.178 ±0.072 ±0.036 
0.323 ±0.094 ±0.019 
0.26 ±0.14 ±0.02 
0.11 ±0.17 ±0.01 
-0.29 ±0.23 ±0.03 


-0.49 ±0.41 ±0.13 
-0.18 ±0.33 ±0.07 

0.10 ± 0.20 ± 0.04 
-0.06 ±0.20 ±0.02 
-0.19 ±0.23 ±0.02 
0.11 ±0:31 ±0.02 

-0.16±g;|| ±0.05 



24 Jean- Philippe Lansberg 



D 
CM 

+ 



(M 



II 



1 

0.5 


-0.5 
-1 



4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 
Pt (GeV) 

Fig. 14. cipronipt for J/'tp fit for \y\ < 0-6- The error bars denote statistical and systematic 
uncertainties added in quadrature. 



ip' polarisation measurement The procedure to obtain aprompt and afromB is 
similar. Weighted simulations of the angular distribution are fit to the data and 
the weight obtained gives the polarisation. The | cos^| distributions in the two ct 
sub-samples are fit simultaneously. As there is no expected radiative decay from 
higher excited charmonia, we are in f act d ealing with direct production. 

Anew, the systematic uncertaintieJSa are much smaller than the statistical un- 
certainties. The values obtained for aprompt and afromB are given in Table. |S| and 
plotted in Fig. [l5l 
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Table 8. Fit results for i/'(25) polarisation, with statistical and systematic uncer- 
tainties (Ref. 



Pt bin (GeV) 


Mean pT (GeV) 


^Prompt 


^ fromB 


5.5 - 7.0 


6.2 


-0.08 ±0.63 ±0.02 


-0.26 ± 1.26 ±0.04 


7.0 - 9.0 


7.9 


0.50 ±0.76 ±0.04 


-1.68 ±0.55 ±0.12 


9.0- 20.0 


11.6 


-0.54 ±0.48 ±0.04 


0.27 ±0.81 ±0.06 



2.5.2. Study of the T{IS) 

The measurements are made in the region ofgr from to 20 GeV and with \y\ < 0.4 
and the data are separated in four pT-hivS^. In Table. El are given the results for 
a and the same values are plotted in Fig. El Our conclusion is that T(15') seems 
to be mostly produced unpolarised. 
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Fig. 16. a for T(1S) fit for \y\ < 0.4. 



Table 9. Fit results for T(1S) po- 
larisation (Ref. i56j. 



Pt bin (GeV) 


a 


0.0 - 3.0 
3.0 - 5.0 
5.0 - 8.0 
8.0 - 20.0 


+0.05 ±0.14 
+0.01 ± 0.14 
+0.10 + 0.17 
-0.12 + 0.22 



2.5.3. New preliminary measurements by CDF for prompt J/tp 

To complete the review of the measurements of quarkonium polarisation done by 
the CDF collaboration, we give in Fig. 1171 the preliminary one for prompt J ftp for 
RUN II with an integrated luminosity of 188 ± 11 pb"^ for J/-0 with 5 < Pr < 30 
GeV and \y\ < 0.6. 
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Fig. 17. Preliminary measurements of the prompt production of J/tp as measured by CDF for 
RUN II (from Ref.CT. 



2.6. PHENIX analysis for J ftp production cross sections 

In this section, we present the first measurements of — > J/4> + X at RHIC 
obtained by the PHENIX experimenJ^ at a cm. energy of — 200 GeV. The 
analysis was carried out by detecting either dielectron or dimuon pairs. The data 
were taken during the run at the end of 2001 and at the beginning of 2002. The data 
amounted to 67 nb~^ for /i^/i^ and'* 82 nb^^ for e+e~. The B-decay feed-down is 
estimated to contribute less than 4% at y^s — 200 GeV and is not studied separately. 
The production is thus assu med her e to be nearly totally prompt, feed-down from 
Xc is expected to exist thougd^^ESI. 

The net yield of Jf-ip within the region 2.80 GeV < M (pair) < 3.40 GeV was 
found to be 46.0 ± 7.4 for electrons, whereas for muons, it was 65.0 ±9.5 J/i/j within 
the region 2.71 GeV < M{pair) < 3.67 GeV. The cross section as a function of p-r 
is shown for the two analyses on Fig. 1181 



"^The difference results from different cuts on the extrapolated vertex position. 
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Fig. 18. The differential cross section measured: left for dielectrons, right for dimuons fRef . JET\ . 



3. Review of contemporary models for production at the Tevatron 
and RHIC 

3.1. Soft- Colour Interaction vs Colour Evaporation Model 

IntroducecPEni as a new way to explain the observation of rapidity gaps in deep 
inelastic scattering at HERA, the Soft-Colour Interaction (SCI) model was also 
applied to quarkonium produc tion , in particular in the context of hadron-hadron 
collisions (Tevatron^^, ^nd LHcj^. 

The main idea of the model is to take into consideration soft interactions occur- 
ring below a small momentum scale, which we shall call here Qq, in addition to those 
considered in the hard interaction through Feynman graphs. Unfortunately, we do 
not have satisfactory tools to deal with such soft QCD interactions. Nevertheless, 
the interesting point of these interaction -emphasised by the SCI model- lies in 
the fact that these matter only for colour, since they cannot affect significantly the 
parton momenta. Therefore, one possibility is to implement them via Monte Carlo 
(MC) event generators by solely exchanging the colour state of two selected softly 
interacting partons with a probability R. From rapidity gap studies, R is close to 
0.5. 



3.1.1. The model in hadroproduction 

To what concerns the hard part, Edin et al. followed a procedure similar to the 
CEM: the prompt cross section to get a given quarkonium Q is obtained from 
the one to get a colour-singlet QQ pair with an invariant mass between 2mQ and 
2mqQ after distributing it between the different states of a family (charmonium or 
bottomonium) . In the CEM, it is supposed -from SU{3) counting- to be one ninth 



28 Jean- Philippe Lansberg 



of the total'* cross section aqQ . 

The latter can be computed by Monte Carlo (MC) event generators, which can 
be coupled to the SCI procedure (possible exchange of colour state between some 
partons). The computation of ctqq can be effectively done either with NLO ma- 
trix ele men ts (which include gluon-fragmentation into a colour-octet pair) or with 
PythitJ^ -including LO matrix elements and parton showers- (which also contains 
gluon-fragmentation into a colour-octet pair) . Choosing Pythia has the advantage 
of introducing even higher contributions which can be significant. To illustrate this, 
we have reproduced a comparison (see Fig. I19|) of two CEM calculations: one with 
NLO matrix-element implementation and another with Pythia Monte Carlo in- 
cluding parton showers, nic was set^ to 1.5 GeV, pj/^p was taken to be 0.5, and 
Pjf,' 0.066. These two choices disagree clearly with simple spin statistics, but are 
necessary to reproduce the normalisation of the data. 
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Fig. 19. Comparison between J/-0 and ip' prompt production cross sections within CEM with 
NLO matrix elements and with Pythia. The data are from CDF and D0 ^'^ (Reprinted figure 
from Ref. IK8l with kind permission of Springer Science and Business Media. Copyright (2002)). 



3.1.2. The results of the model 

The cross section to produce the colour-singlet heavy-quark pair (see Eq. Q)) is com- 
puted here with Pythia which is coupled to the SCI model: once the selected par- 
tons -with a probability R- have exchanged their colour, only the events producing 
a colour-singlet heavy-quark pair are retained. This directly gives the singlet-state 
cross section, whereas in the CEM ■^(^qq is used. It is also integrated in the mass 

^The epithet total refers to the colour: singlet plus octet configurations, not to a possible integration 
over Pj-. The dependence on Pj- of a is implied and not written down to simplify notation, 
''these are the values to be taken to reproduc J^^ l fixed-target measurements. 
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region between 2mQ and 2mqQ. The cross section to produce one given quarkonium 
Q is obtained using: 

'^Q = -p CToniuTm (19) 

with = where Ji is the total angular momentum and rii the main quantum 

number. pQ is thus effectively replaced by and is not free to vary anymore. 

The effect of SCI can be either to turn a colour-singlet pair into a colour-octet 
one or the contrary. The latter case is important, since it opens the possibility of 
gluon fragmentation into a quarkonium at order (see Fig. I20|) '^. This explains 
why the model reaches a reasonable agreement with data to what concerns the 
slope vs pt - the same is equally true for the CEM. Note that the effect of SCFs 
on the cross section depends on the partonic state (through the number of possible 
SCI's), and thus on the transverse momentum of the produced quarkonium. The 
final colour-singlet pT slope is slightly steeper than the initial QQ one and in better 
agreement with data (compare Fig. [T^ and Fig. |2] (left)) 



SCI 




Fig. 20. Gluon fragmentation process (part of NLO contributions) giving a singlet state through 



The very interesting point of this approach is its simplicity. R is the only free 
parameter and is kept at 0.5, which is the value chosen to reproduce the rate of 
rapidity gaps. However, the cross section depends little on it. Putting R to 0.1 
instead of 0.5 decreases the cross section by 30%. 

On the other hand, in Eq. H19() , the ^ suppression for radially excited states does 
not fundamentally follow from spin statistics; it is solely motivated by the ratio of 
the leptonic decay width. 

'^The same thing happens implicitly in the CEM since one ninth of ctqq gives colour-singlet pair 
irrespectively of the kinematics and thus of which diagram in the hard part contributes the most. 
This can also be compared to the Colour-Octet Mechanism as one can extract a fragmentation 
function from the simulatioil^^. 
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Another drawback is the strong dependence upon the heavy-quark mass*^. This 
is mainly due to modified boundary values in the integral from 2rnQ to 2mqQ , much 
less to a change in the hard scale and in the mass value entering the matrix elements. 
Changing rric from 1.35 GeV to 1.6 GeV decreases the cross section by a factor of 3, 
and from 1.35 GeV to 1.15 GeV increases it by a factor of 2. The other dependences 
(on Aqcd, PDF sets, ...) are not significant. 
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Fig. 21. SCI calculations for prompt production of J/tp, i>' (left) and T(1S) and T(25) (right). 
The dashed lines correspond to the LO matrix-element alone, the solid one to the LO matrix- 
elements with parton showers (Reprinted figures from Ref. ltibl with permission of American Physical 
Society. Copyright (1997)). 



3.1.3. Improving the mapping 

As one can imagine, physically, during the transition between the QQ pair - 
produced by the hard interaction with an invariant mass between 2toq and 2mqQ- 
and the physical quarkonium Q, this invariant mass is likely to be modified. It is 
then reasonable to suppose that it is smeared to a mass m around its initial value 
mQQ with such distribution: 

/ {niQQ — m) \ 
Gsnie{mQQ,m) = exp I I . (20) 

\ sine / 

Supposing that the width of the quarkonium resonances tends to zero, the prob- 
ability that a pair with an invariant mass toqq gives a quarkonium Qi is given by: 

. (21) 

stih with r, = 

This smearing enables quark pairs with invariant masses above the heavy-meson 
threshold 2mqQ to produce stable quarkonia, but the reverse is true also, with a 
probability A{mQQ) equal to ^erfc ( ^"^J^J^^^'^ 



"^The same dependence exists for the CEM. 
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The cross section to get a quarkonium Qi is thus given Ref. by 



^ . da, 
dm 



2m. 



'^'^QQTi 7(1 - -4(TOQQ))P,(mQQ), (22) 



Q 



with calculated either within CEM (with \ from S'C/(3)) or within SCI. 

Using this procedure, the production ratio ip' over J/ip is better reproduced^ 
than for simple spin stat istic s as well as the different components of the prompt 
J/il) as measured by CDeESI On the other hand, due to the proximity in mass of 
the Xci and Xc2, their production ratio is not affected and remains 3 : 5 slightly in 
contradiction with the CDF measurement consistent with a ratio of 1. Finally, 
this refinement affects the cross sections only by a factor of 20%, what is well within 
the model uncertainty. 



3.2. NRQCD: including the colour-octet mechanism 

In 1992, Bodwin et al. considered new Fock-state contributions in order to cancel 
exactly the IR divergence in the light-hadron decays of Xci (and he) at LO. This 
decay proceeds via two gluons, one real and one off-shell, which gives the light- 
quark pair forming the hadron. Whe n the first glu on becomes soft, the decay width 
diverges. The conventional treatmenlP3^3^ESI, which amounts to regulating the 
divergence by an infrared cut-off identified with the binding energy of the bound 
state, was not satisfactory: it supposed a logarithmic dependence of ^''(0) upon 
the binding energy. They looked at this divergence as being a clear sign that the 
factorisation was violated in the CSM. 

Their new Fock states for Xc were e.g. a gluon plus a cc pair, in a '^Si configura- 
tion and in a colour-octet state. The decay of this Fock state occurred through the 
transition of the coloured pair into a gluon (plus the other gluon already present in 
the Fock state as a spectator). This involved a new phenomenological parameter, Hg, 
which was related to the probability that the cc pair of the Xc be in a colour-octet 
S'-wave state. The key point of this procedure was that a logarithmic dependence 
on a new scale A - a typical momentum scale for the light quark - appeared nat- 
urally wit hin th e effective field theory Non-Relativistic Quantum Chromodynamic 

(nrqcd|SIZZI. 

This effective theory is based on a systematic expansion in both as and v, 
which is the quark velocity within the bound state. For charmonium, ~ 0.23 and 
for bottomonium ~ 0.08. One of the main novel features of this theory is the 
introduction of dynamical gluons in the Fock-state decomposition of the physical 
quarkonium states. In the case of S'-wave orthoquarkonia {^Si), we schematically 
have, in the Coulomb gauge: 

\QQ)^0{l)\QQ['s['^])+Oiv)\QQ['pfg])+0{v^)\QQ['S^,'^g]) 
+ 0{v^)\QQ[ 'Si'^'^gg]) + 0{v')\QQ[ ^D^^^'^gg]) + ... 
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whereas for P-wave orthoquarkonia {'^Pj), the decomposition is as follows 

\Qq,) ^0{l)\QQ[^P^j'^]) + 0{v)\QQ[^s['^g]) + . . . (24) 

In these two formulae, the superscripts (1) and (8) indicate the colour state of the 
QQ pair. The 0{v^) factors indicate the order in the velocity expansion at which 
the corresponding Fock state participates to the creation or the annihilation of 
quarkonia. This follows from the velocity scaling rules of NRQCD (see e.g. Ref . I77|l . 

In this formalism, it is thus possible to demonstrate, in the limit of large quark 
mass, the factorisation between the short-distance - and perturbative - contri- 
butions and the hadronisation of the QQ, described by non-perturbative matrix 
elements defined within NRQCD. For instance, the differential cross section for the 
production of a quarkonium Q associated with some other hadrons X reads 

da{Q + X) =Y,dHQQ?'''''L[j''^] + X){0^[''^+'L'^j''^]), (25) 

where the sum stands for S, L, J and the colour. 

The long-distance matrix element (LOME) {O^ takes account of the 

transition between the QQ pair and the final physical state Q. Its power scaling rule 
comes both from the suppression in v of the Fock-state component '^^] in 

the wave function of Q and from the scaling of the NRQCD interaction responsible 
for the transition. 

Usually, one defines 0^[^^^^l[^'^^] as the production operator that creates and 
annihilates a point-like QQ pair in the specified state. This has the following general 
expression 

[25+1^(13)] ^ ^t^v, f ^ 5: IQ + X)(Q + Xl") ^^Kx 

\ X J, J (26) 

where® ip and x ^re Pauli spinors and the matrix K is a product of colour, spin 
and covariant derivative factors. These factors can be obtained from the NRQCD 
lagrangiarf^. For instance, ^s[^^] = X^''^V'('2gfls)'0^<?X- However some transi- 
tions require the presence of chromomagnetic (AL = and AS* = ±1) and chro- 
moelectric (AL = ±1 and AS — 0) terms for which the expressions of K are more 
complicated. 

On the other hand, the general scaling rules relative to the LDME's^^ give: 

(©2 [2-5+1 L'-j^^^]) y3+2L+2E+4M^ (27) 

where E and M are the minimum number of chromoelectric and chromomagnetic 
transitions for the QQ to go from the state [^^^^L^j'^''] to the dominant quarkonium 
Q Fock state. 



"The second line of Eq. 1261 is nothing but a short way of expressing this operator. 
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The idea of combining fragmentation as the main source of production with 
allowed transitio ns b etween a Xc to a ^S'l in a colour-octet state, was applied by 
Braaten and YuarP^. Indeed, similar formulae to the one written for fragmentation 
within the CSM can be written for fragmentation functions in 

I?,^s(z,m) =5]dp,+,^a.s,](z,/i)(O^P^+'4'''^])> (28) 

where accounts for short-distance contributions and does not depend on 

which Q is involved. But since the theoretical predictions for prompt J /^p production 
did not disagree dramatically with data and since there was no possible Xc decay 
to a possible enhancement of the Xc cross section by colour-octet mechanism 
(COM) was not seen at that time as a key-point both for J /-ip and ij}' production. 

However, in the case of J /"ip and -0' production, COM could still matter, but in 
a different manner: fragmentation of a gluon into a ^-Pj*"* is possible with solely one 
gluon emission and fragmentation into a ^sf'^ requires no further gluon emission 
(at least in the hard process -described by d^.^{z,^)- and not in the soft process 
associated with O). Concerning the latter process, as two chromoelectric transitions 
are required for the transition \ccgg) to |cc), the associated LDME {O"^^ S^^'']) was 
expected to scale as rn^v^ . In fact, d,^ ^s),, the contribution to the fragmentation 

function of the short- dist ance process g S^^^ was already known since the paper 
of Braaten and Yuarl^ and could be used here, as the hard part (i[.](z,/j,) of the 
fragmentation process is independent of the quarkonium. 




p, (GeV) 



Fig. 22. Differential cross section versus pT of the CSM (dotted: LO; dashed: fragmentation and 
LO) production and of the COM fragmentati on ( solid curve) to be compared with CDF preliminary 
measurements of the direct production of i/j'tEM (Reprinted figure from Ref. 1781 with permission of 
American Physical Society. Copyright (1995)). 

In a key paper, Braaten and Fleming combined everything together to cal- 
culate, for the Tevatron, the fragmentation rate of a gluon into an octet '^Si that 
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subsequently evolves into a ip' . They obtained, with^ (O'^^^fU^^f' ]) = 4.2 x lO'^ 
GeV'^, a perfect agreement with the CDF preliminary datasi^ compare these pre- 
dictions in Fig. 1221 with the CSM fragmentation ones in Fig. [S] (right). 

Following these studies, a complete s urvey on the colour-octet mechanism was 
made in two papers by Cho and Leibovic n' ^ 1 "^ The main achievements of these pa- 
pers were the calculation of colour-octet P-state contributions to V', the predictions 
for prompt and direct colour-octet J/tp production with also the Xc feed-down cal- 
culation, the first predictions for T and the complete set of 2 — > 3 parton processes 
like ij —fcc+k^ all this in agreement with the data. 

3.2.1. Determination of the LDME's 

Here we expose the results relative to the determination of the Long Distance Matrix 
Elements of NRQCD necessary to describe the production of "^Si quarkonium. We 
can distinct between two classes of matrix elements : the colour singlet ones, which 
are fixed as we shall see, and the colour octet ones which are fit to reproduce the 
observed cross sections as a function of px- 

In fact, as already mentioned, NRQCD predicts that there is an infinite number 
of Fock-state contributions to the production of a quarkonium Q and thus an infinite 
number of LDME. Practically, one is driven to truncate the series; this is quite 
natural in fact since most of the contributions should be suppressed by factor of at 
least C(^), where v is the quark velocity in the bound state. 

For definiteness, the latest studies accommodating the production rate of "^^i 
states retain only the colour singlet state with the same quantum numbers as the 
bound state and colour octet P-wave states and singlet iS-wave. In this context, the 
CSM can be thought as a further approximation to the NRQCD formalism, where 
we keep solely the leading terms in v. 

Colour-singlet LDME In this formalism, factorisation tells us that each contribu- 
tion is the product of a perturbative part and a non-perturbative matrix element, 
giving, roughly speaking, the probability that the quark pair perturbatively pro- 
duced will evolve into the considered physical bound state. If one transposes this 
to the CSM, this means the wave function at the origin corresponds to this non- 
perturbative element. This seems reasonable since the wave function squared is also 
a probability. 

Yet, one has to be cautious if one links production processes with decay pro- 
cesses. In NRQCD, two different matrix elements are defined for the "colour singlet" 
production and decay, and they are likely to be different and independent. 

The only path left to recover the CSM is the use of a further approximation, the 
vacuum saturation approximation. The latter tells us how the matrix element for 

'This corresponds to a suppression of 25 compared to the "colour-singlet matrix element" 
(^Qi^(2S) [35(1)])^ which scales as m^v^. This is thus in reasonable agreement with a ti^ suppression. 
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the decay is linked to the one for the production. This enables us to relate the wave 
function at the origin appearing in Tu (Eq. (3.1) of Ref. ED to the colour-singlet 
NRQCD matrix element for production. This gives: 

{O'^' [3^(1)]) ^ 18 1^(0)1' + 0{v^). (29) 

The conclusion that could be drawn within NRQCD is that the extraction of 
non-perturbative input for production from the one for decay is polluted by factors 
of O(u^), this is also true for extraction from potential models. 

In Table. 1101 we give the colour-singlet LDME for the J /ip and the -0'. The 
result for the different potentials are deduced from the solutions of Ref. |S2J These 
LDME's -up to a factor 18- are those that are to be used in CSM calculations. The 
values differ from the one used in Refs. [Tfl I28L 1291 because of modifications in the 
measured values of r«, NLO QCD corrections to Tu and also in the potential used 
to obtain the wave function at the origin. 



Table 10. Colour-singlet LDME for the J/i/> and the j/i' determined from the leptonic decay width 
and from various potentials. The error of 0(v*') as shown in Eq. I29i should be implied. Values are 
given in units GeV'^. 



Q 


Leptonic decay 


BT Potential Ref. 83 84 


POW potential Ref. 85 


Cornell Potential Ref. 86 8^^ 




l.lit 0.1 


1.16 


1.4 


2.06 


^' 


0.7± 0.1 


0.76 


0.76 


1.29 



In Table. im we expose the results of Ref. ISHlconcerning the colour-singlet LDME 
for T(nS'), i.e. (O^t"^' pS*!^^]). The error quoted for the value from potential models 
expresses the variation of the latter when passing from one to another. 

Table 11. Colour-singlet LDME for the 
T(nS) determined from the leptonic de- 
cay width and from potentials models. 
Values are given in units GeV^ (Ref. IS^ . 



Q 


Leptonic decay 


Potential Models 


T(1S) 
T(2S) 
T(3S) 


10.9 ± 1.6 
4.5 ±0.7 
4.3 ±0.9 


10.8 ± 5.5 
5.0 ± 1.8 
3.7± 1.5 



Colour-octet LDME's As said above, three intermediate colour-octet states are 
currently considered in the description of ■^Si production. These are ^Sq^\ '^Pq^^ 
and ■^S'j;^''. The corresponding LDME's giving the probability of transition between 
these states and the physical colour-singlet ^Si state are not known and are to be 
fit to the data. 

Unfortunately, the perturbative amplitudes to produce a ^S^^^ or ^P^^^ have 
the same px slope and their coefficient cannot be determined apart. Therefore, 
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one defines k as the ratio between these two amplitudes. From it, one defines the 
following combination 

M^CSI,'\'P^'^) ^ {0^['sl,'^])+k^-^^^^, (30) 

which is fit to the data. In the following, we expose the results obtained by different 
analyses using various PDF set and parameter values. In the following tables, the 
first error quoted is statistical, the second error, when present, reflects the variation 
of the fit LDME when the renormalisation and factorisation scales is set to ^ — 
1/2 y^py + Ami and to 2 A/p|r+4TOf. The agreement with the data being actually 
good, there is no real interest to plot the cross sections given by the fits. 



Table 12. Fit values of J/ip production LDME's from at the Tevatron. Values are given 
in units lO^^ GeV^ (Numbers from Ref.l5?7t. 



Reference 


PDF 






k 


Cho-Leibovich Ref. 80 


MRS (DO) Ref. 89 


0.66 ±0.21 


6.6 ± 1.5 


3 


Beneke-Kramer Ref. 1901 


CTEQ4L Ref. 91 
GRV-LO(94) Ref. 92 
MRS(R2) Ref.;i3. 


1.06 ± 0.14+Q^i^ 
1.40±0.221q;^^ 


4.38 ± 1.15+o;74 
3.90± l.uj;^;;^^ 
10.9± 2.07l^;2e 


3.5 


Braaten-Kniehl-Lee Ref. 


MRST-LO(98) Ref. 155 
CTEQ5L Ref. ESI 


0.44 ± 0.07 
0.39 ±0.07 


8.7 ± 0.9 
6.6 ± 0.7 


3.4 


Kramer Ref. 50 


CTEQ5L Rei.W 


1.19 ±0.14 


4.54 ± 1.11 


3.5 



Table 13. Same as Table [T^ for ip(2S) production. Values are given in units 10 ^ GeV'^ (Numbers 
from Ref. EOJ 



Reference 


PDF 




Mf"')(i5(f',aPo'«') 


k 


Cho-Leibovich Ref. ESI 


MRS (DO) Ref. iM 


0.46 ±0.1 


1.77 ± 0.57 


3 


Beneke-Kramer Ref. 1901 


CTEQ4L Ref. TO 
GRV-LO(94) Ref.ini 
MRS(R2) Ref.ini 


0.44 ± 0.08+o;24 
0.46 ± 0.08+023 
0.56±0.1ltJ!;g 


1.80± 0.56+030 
1.60±0.5ll»to 
4.36±0.96ti;iJ 


3.5 


Braaten-Kniehl-Lee Ref. |94 


MRST-LO(98) Ref. gF 
CTEQ5L Ref. 96 


0.42 ±0.1 
0.37 ± 0.09 


1.3 ± 0.5 
0.78 ±0.36 


3.4 


Kramer Ref. 50 


CTEQ5L Ref. J£_ 


0.50 ± 0.06 


1.89 ± 1.11 


3.5 



In the following tables (Tables 12-16), the results of Cho and Leibovich are for 
Pt > 3.5 GeV on the data of^^ and the ones of Braaten et al. are for px > 8 GeV 
and rub = 4.77 GeV was chosen. 



Table 14. Same as Table ini for T(15) production. Values are given m units 10-2 GeV^. 



Reference 


PDF 


(0^(1S)[35(8)]) 


A^T(15)(l5(8,^3p^(8)) 


k 


Cho-Leibovich Ref. .80 


MRS(DO) Ref. 89. 


0.6 ±0.2 


0.4 ± 0.5 


5 


Braaten Fleming Lei- 
bovich Ref. 1551 


MRST-LO(98) Ref.|5^ 
CTEQ5L Ref. El 


o.4±o.7;i.- 

2.0±4.i;0j 


20.2 ± 7.8+g'^g'' 
13.6±6.8l!,°g** 


5 
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Table 15. Same as Table [T^ for T(2S') production. Values are given in units 10 ^ GeV''. 



Reference 


PDF 






k 


Cho-Leibovich Ref. 80 


MRS(DO) Ref. 89 


0.4 ± 0.1 


0.5 ±0.4 


5 


Braaten Fleming Lci- 
bovich Ref.W 


MRST-LO(98) Ref. 12^ 
CTEQ5L Ref.l96l 


17.4 ±6.4+^;^ 
16.4± 5.7tg } 


-9.5±ll.l+^;» 
-W.%±9.7lli 


5 



Table 16. Same as Table IT2I for T(3S') production. Values are given m units 10^2 GeV^. 



Reference 


PDF 


(ClT(3S)[35t«)]^ 




k 


Braaten Fleming Lei- 
bovich Ref. 1551 


MRST-LO{98) Ref. 95 
CTEQ5L Ref. .96. 


3.7± 2.1+i;;^ 
3.6 ± 1.9j;i;3 


7.5±4.9+^;| 
5.4±4.3j:^;^ 


5 



3.2.2. Polarisation predictions 

A straightforward and unavoidable consequence of the NRQCD sohition to the ij]' 
anomaly was early raised by Cho and WisJ^ the il)' , produced by a fragment- 
ing (and real) gluon through a colour octet state, is 100% transversally polarised. 
They in turn suggested a test of this prediction, i. e. the measurement of the lepton 
angular distribution in ijj' — > £^i^, which should behave as j^^(V'' ^ oc 
(1 + acos^d), with a=l for 100% transversally polarised particles since the spin 
symmetry of NRQCD prevents soft-gluon emissions to flip the spin of QQ states. 

In parallel to the extraction studies of LDME, the evolution of this polarisation 
variable a as a function of pr was thus predicted by different groups and compared 
to measurement of the CDF collaboration (see section . 

If we restrict ourselves to the high-pr region, where the fragmenting gluon 
is transversally polarised, the polarisation can only be affected^ by v corrections 
(linked to the breaking of the NRQCD spin symmetry) or as corrections different 
than the ones already included in the Altarelli-Parisi evolution of the fragmentation 
function D. Indeed, emissions of hard*^ gluons are likely to fl ip t he spin of the QQ 
pair. These corrections have been considered by Beneke et a/!^. 

In Fig. [531 we show the various polarisation calculations from NRQCD for 
prompt J/ip, direct ip' and T(IS') with feed-down from higher bottomonium taken 
into account. The least that we may say is that NRQCD through gluon fragmenta- 
tion is not able to describe the present data on polarisation, especially if the trend 
to have a < at high pT is confirmed by future measurements. 

Motivated by an apparent discrepancy in th e hierarchy of the LDME's for J/ip 
and -ip' J Fleming, Rothstein and Leibovicl^^ proposed different scaling rules be- 
longing to NRQCDc. Their prediction for the cross-section was equally good and 
for polarisation they predicted that a be close to 1/3 at large pr- 



^We mean a 7^ +1. 

''with momenta higher than A^nrqcd- 
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inclusive T(1S) 



20 25 30 
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Fig. 23. NRQCD calculation of a{pT) for (left) J/ip (Braaten et al: Ref.^U, (middle) i/;' (Lei- 
bovich: Ref.|HHl Braaten et al: Ref.|Hll Beneke et al: Ref.EJ and (right) T(15) (Braaten-Lee: 
Ref.|nni-(left reprinted figure from Ref.lMlwith permission of American Physical Society. Copyright 
(2000); from Ref.lHl middle and right from Ref.lTTt 



This latter proposal nevertheless raises some questions since the very utility 
of the scaling rules was to provide us with the evolution of the unknown matrix 
elements of NRQCD when the quark velocity v changes, equally when one goes 
from one quarkonium family to another. Indeed, a LMDE scaling as 0{v'^) may be 
bigger than another scaling as 0{v^) since we have no control on the coefficient 
multiplying the v dependence. On the other hand, a comparison of the same LDME 
for a charmonium and the corresponding bottomonium is licit. Now, if the counting 
rules are modified between charmonia and bottomonia, the enhancement of the 
predictive power due to this scaling rules is likely to be reduced to saying that the 
unknown coefficient should not be that large and an operator scaling as 0{v^) is 
conceivably suppressed to one scaling as 0{v^), not more than a supposition then. 



3.3. kx factorisation and BFKL vertex 

If one considers the production of charmonium or bottomonium at hadron colliders 
such that the Tevatron, for reasonable values of pr and of the rapidity y, it can 
be initiated by partons with momentum as low as a few percent of that of the 
colliding hadrons. In other words, we are dealing with processes in the low Bjorken 
X region. In that region we usually deal with the BFKL equation, which arises from 
the resummation of as\og{l/x) factors in the partonic distributions. This process 
of resummation involves what we can call Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
effective vertices and the latter can be used in other processes than the evolution 
of parton distributions. 

On the other hand, the kr factorisation approadP^ ^ I ^'^^^ l ^'^^'^ l ^'^^^ l ^^^^ l ^^^^ which 
generalises the coUinear approach to nonvanishing transverse momenta for the initial 
partons, can be coupled to the NLLA vertex to describe production processes, such 
as those of heavy quark or even quarkonium. 

This combination of the fc^ factorisation for the initial partons and the NLLA 
BFKL effective vertex -'^^^ for the hard part can be thought as the natural frame- 
work to deal with low x processes since its approximations are especially valid in 
this kinematical region. As an example, very large contributions from NLO in the 
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coUinear factorisation are already included in the LO contributions of this approach. 
A typical case is the fragmentation processes at large px in quarkonium production. 

Since this approach is thought to be valid for low x processes, but still at a 
partonic scale above Aqcd , this can be also referred to as the semi-hard approach. 
One can find an useful review about the approach, its applicat ions and its open 
questions in the two following papers of the small-x 

coUaboratioiPHCni. 



3.3.1. Differences with the collinear approach 

Practically, compared to the collinear approach, we can highlight two main differ- 
ences. First, instead of appealing to collinear parton distributions, we are going to 
employ the unintegrated PDF J-{x^ qj^) for which there exist different parameter- 
isations, exactly as for the collinear case. The discussion of the difference between 
these is beyond the scope of this review. 

These are however related to the usual PDF's by: 

xg{x,p?)^ r M:^(x,qT). (31) 
Jo <3t 

Secondly, the hard part of the process is computed thanks to effective vertices 
derived in Rcf. 107 or following the usual Feynman rules of pQCD using an extended 
set of diagram -'^-'^^ due to the off-shell ness of the i-channel gluons and then using 
what is called the Collins-Ellis "trick'^^^ which consists in the following replace- 
ment: e''(qi)e*''((7i) — j^^^- Let us present here the first approach and give the 

expression for the QQ production vertex 

= -gi52«(fci) (T'^T' B{qi,q2, fci, fe) - TT^B^Cgi, 92, fe, fci)) v{k2), (32) 

where a, b are the gluon colour indices, T"^ — ^A^, fci, k2 the (on-shell) heavy- 
quark and antiquark momenta, qi, q2, the (off-shell) gluon momentum, and gi, g2 
the strong coupling constant which is evaluated at two different scales, q^ and 
respectively. The expression for B{qi, q2, fci, fe) (and i3^) is a sum of two terms 



„t , , \ _^i_L-^i_L-TO + 7/3?+ ^(g2,gi) 

e(<Zl,g2,fcl,fc2)=7 (^^_fc^)2_^2 7 - (fc^+fc,)2 ^ 

S-(,i,,„fc„fcO =7+iii^%±^7- ^'^r'^'^^-'^^) 



[qi - k2Y - rn^ (fci -I- k2Y 



(33) 



where T+- Hq2.qi) = 2{q^- q2f ~ 2q+n-P + 2q^n+f' - 2ql + 2q^ 

The first corresponds to the contribution of B in Fig. |^ the second, with an 
^ propagator, is linked to a transition between two (off-shell) i-channel gluons 
(Reggeons) and a (off-shell) gluon which subsequently splits into the heavy quark 
(C in Fig. EH). 

This term is not only derived from a triple gluon vertex (D of Fig. 124(1 . The 
complication is expected since the gluons are off-shell. Indeed if we want to deal 
with on-shell particles to impose current conservation, we have to go back to the 
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Fig. 24. Different components in the QQ production BFKL vertex. 

particles -here the initial hadrons- which have emitted these i-channel gluons. Those 
can directly emit the third gluon in addition to the i-channel gluon (E of Fier. I24|l . 
similarly to a Brcmfitrahlmig contribution. This kind of emission gives birth to the 
two terms: — 2g? — - + 2qi — -r. The factors g? and g? account for the fact that 
there is only one t-channel gluon in this type of process, the denominator comes 
from the propagator of the particle which has emitted the s-channel gluon. 

The vertex of Eq. H33|) has been used successfully for the description of open 
beauty production-'^ "'^-'^ and was also shown to give a large contribution to the pro- 
duction of Xci in the Colour-Singlet approacti I 

In the case of tp and T production in the Colour-Singlet approach, we need to 
consider the vertex g*g* — > QQg rather than g*g* QQ to conserve C-parity. This 
would correspond to NNLLA corrections and these are not known yet. However, 
the reason of this complication (C-parity) is also from where one finds the solution 
since it sets the contributions of the unknown diagrams to zero when projected on a 
colour-singlet ^Si QQ state. The expressions are finally the same with the addition 
of one gluon emission on the quark lines (see R.ef. IllHj l. 

The cross section is then obtained after the integration on the transverse mo- 
mentum of the gluons' and on the final state momenta [P is the quarkonium mo- 
mentum, k that of the final state gluon). As in the usual CSM, the heavy-quark 
pair is projected on a colour-singlet state and their relative momentum is set 



'The fractions of momentum carried by the gluons xi_2 = are automatically integrated during 

^1,2 

this integration. 
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to zero (fci — k2 and P = ki + . This givesj 
1 1 fd^Pd^k 



1 If d P d k -> - 

<^PiP2^QgX = lg(27r)4 8^ J "^"^'^^9iTd^'f2T'^^('fiT - g2T - kiT - fc2T)x 



■F(xi,giT) ■F(xi,g2T) I ■ .,2 

(34) 

where - a and b being the colour indices of the t-channel gluons, i and j those of the 
heavy quarks ((3i,3j|l) = ^) and TZ{0) the radial part of the Schrodinger wave 
function at the origin (in position space) - 



■^"'^E E ^(0,0;l,5,|l, J,)(3z,3j|i; 



7^(0) 

■X 



X Tr 



Different projections can be used if one studies other quarkonia than the '^5'i. 
This formalism can be also combined with the COM by projecting the heavy-quark 
pair on a colour-octet state and introducing a LDME to give the probability for the 
non-perturbative transition into the physical quarkonium, exactly as in the collinear 
approach presented in section 13.21 Of course, as we shall see in the results, these 
LDME's will have modified values compared to the ones of the collinear fits. 



3.3.2. Results for J / ^ , ip' and T(IS) 

Let us first present the result of Hagler et al^^^ who have used the KMS unin- 
tegrated PDF-H^. Since the hard part of the process is modified compared to the 
collinear case, it is not surprising that the ratio k of the cross section of ^Sq^^ and 
^Pq^^ is different, the slopes are still similar, though; an independent fit would not 
again distinguish the two LDME's associated with these processes. The ratio k lies 
between 6 at low Pt and 4.5 and high Pt- They did the fit with these two contribu- 
tions alone; their results are given in Table. El The first combination (only ^S^^^) 
is plotted in Fig. Esl 



Table 17. Same as Table. IT2I for J/^p production in the kx factorisation ap- 
proach and the NLLA BFKL vertex. Values are given in units 10~^ GeV"^. 





PDF 






k 




Hagler et al. 1 


KMS Ref. 114 0.032 ± 0.012 


7.0 ±0.5 


5 


only^5g 


Hagler et al. 2 


KMS Ref. 114, 0.05 ± 0.012 


6 ±0.5 


5 


only^Pj 



^To make connection with the "trick" of Collins and EUis^" as presented i note the presence 

in the denominators of q^j, and which would come from the following replacement of the t- 

channel gluon polarisation vectors: ei^(qi)€*'^(qi) = -[^^p|3-- The vectors in the numerators 

would appear e.g. in the three first terms of F"' 
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Fig. 25. J 1^ production cross sections at ^fs = 1.8 TeV for the CSM in the collinear approach 
(LO + fragmentation), in the fcy factorisation approach for the Colour-Singlet contribution and 
the Colour-Octet (^S^*' and ^Sq*') compared to the direct measurement of CDF '^Si (Reprinted 
figure from Ref. 11131 with permission of American Physical Society. Copyright (2001)). 



To what concerns the colour-singlet contribution, as can be seen on Fig. 1531 it 
is more than one order of magnitude larger than the LO CSM contribution, and 
also larger than fragmentation-CSM contribution. The same trend is confirmed by 
the '0' cas This can be partially explained (up to a factor 2.5) by a genuine 
different choice for the scale at which should be evaluated (cf. the factors g\ and 
52 in Eq. 

On the other hand, the colour-octet LDME (O'^/'^pS'f ^]) is thirty times 
smaller than in the collinear fits (compare Table. ^1 and Table. I17|l whereas 
M^.^ ^ Q- s'^^ ^'^ p'^^ ) is similar. In Ref. 11131 it is however emphasised that 
(O^/^psf)]) =0 would give a much worse fit, whereas, according to BaranoA^^, 
{p^l'^fsf^\, can be set'' to with the unintegrated PDF of Ref . lUJj The latter 
analysis of Baranov is dominated by the Colour-Singl et co ntribution and due to the 
longitudinal polarisation of the initial off-shell gluonJH^, the polarisation parame- 
ter a should be negative and close to -1 as soon as Pt reaches 6 GeV both for J/?/; 
and V''- The same results for J/-0 and "0' as presented in Fig. 3 of Ref. lllSl are irrel- 
evant since they do not take into account the weights of the different contributions 
in the particle yield. The polarisation of the quarkonia from isolated colour-octet 
channels is not measurable and it is obvious from the cross-section plots (Fig. 1 and 
Fig. 2 of Ref. 118) that none of the colour-octet channels can reproduce the cross 
section alone, they should be combined. 

To conclude this section, let us mention the combination of the factorisation 



''In fact, the values obtained in this work are not fr om fit. The theoretical uncertainties highlighted 
there were too large to make a fit mcaningfuW^. 
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approach with the colour-octet gluon fragmentation by Saleev et al\ \ The hard 
part which is considered is g*g* g where the g* are reggeised gluons (again 
distributed according to an unintegrated PDF) and where the gluon, g, subsequently 
fragments into a quarkonium via the colour-octet state ^S'J^\ exactly like in the 
collinear approach. The LDME values used were the one from the collinear fij^ 
(see Table. 1121 : (O^'/'^psf ^]) = 0.44 x IQ-^ GeV^ and (C^'pS'f ]) = 0.42 x lO'^ 
GeV'^. The agreement is reasonable co ntras ting with the conclusion from the LO 
kx factorisation analysis of Hagler et a/!^^^' which requires a strong suppression of 
the ^Sf ^ channel. 



3.4. Durham group: enhanced NNLO contributions 

As already said, ^Si quarkonia produced by gluon fusion are necessarily accompa- 
nied by a third gluon in the final state. Indeed, it is required for C-parity conserva- 
tion and in the case of semi-inclusive reaction, as the ones we have been considering 
so far, this gluon cannot come from the initial states (see Fig. ISEl (left)). 

As we have seen, the classical description -through the CSM- of this kind of 
production (especially at LO) in QCD severely underestimates the p roduction rates 
as measured at the Tevatron and even at RHIC. In their wor ll ^ 1 V.A. Khoze et 
al. considered the special case where this third gluon, attached to the heavy-quark 
loop, couples to another parton of the colliding hadrons and produces the gluon 
needed in the final state. Indeed, it is likely that the large number of possible 
graphs -due to the large number of available gluons at large s- may compensate 
the as suppression. The parton multiplicity n behaves like logs and this process 
(see Fig. 1211 (right)) can be considered as the LO amplitude in the BFKL approach 
whereas it is NNLO in pQCD. 



3.4.1. Integrated cross section 

Since the two t-channel gluon off the quark loop are now in a colour-octet symmetric 
state, the real part of the amplitude is expected to dominated by its imaginary part 
(in Fig. 1261 (right) the two quarks and the gluon is the s-channel are then put on- 
shell as well as the remaining quark entering the Q). One then chooses to work in 
the region where the rapidity between the Q and the fi nal gluon p is large {i.e. when 
s ^ Mq) since it should dominate and one get the following expressions for 
the imaginary part of the amplitude for the two possible Feynman graphs (the two 
other ones with the loop momentum reverted give a factor two): 



9mA" =-d"''"5(47ra,)5/2x 
8 



Tr(^(pi)(f + mQ)mM^ / + ^nQ)M£ - A + ^q)) ^^^^ 



2^,s((£ - p,)^ A2)((g -I- ly Xl) 
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(a) (b) 




Fig. 26. (left) Usual LO pQCD in the CSM; (right) NNLO pQCD or LO BFKL contributions. 
In both cases, the gg — > Qg sub-process is shown in bold and the two quarks entering the Q are 
on-shell. 



2 Tr( ^(pi)(f - £ + mQ)M^ + mQmP)Mi - ^ - 7^1 + mg)) 
2ns{{^ - £ - p,Y - A^)((g + £? A^) 

(37) 



5 is a quantity related, up to some known factors, to |'0(O)p through the leptonic 
decay width (g^ = f "f a ) g^^h that the '^Sx vertex reads f (f + MQ)e*{P). \„ 
is an effective gluon mass to avoid logarithmic infrared singularity, motivated by 
possible confinement effects. 

The partonic differential cross section thus reads: 

- ' ' (38) 



dP| f 67rs2 ' 

with A — 2(j4° + A''). The hadronic cross section is obtained with the help of 
da f dx2 , , , .da{s,P^) 



dydPf 



/ xig{xi)x2gix2) — — , 39 

J X2 dPf 



s — SX1X2 with the collision energy in the hadronic frame, y is the rapidity of 
the Q also in the latter frame, g{xi) are the gluon PDF. 

Due to the low-x behaviour of the PDF, the main contribution to the integral 

comes from the lowest value of X2, that is X2 — — — — -e^^ (and thus xi ~ 
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— — — -e^). For reasons exposed in Ref. I120l the a;2-integration' region is to be 
extended over the whole kinematically available rapidity interval Ay, this gives 

= .i(y)5(.i(y)).,(y)g(x.(y))Ay^^il%0ll^, (40) 



with xi.2{y) = V^'^Q_+^? e±y and 



2 

T S 

■^■m /it'-' 



with Xmaa; Set to 0.3 to exclude contributions when the third gluon couples to 
partons with a; > 0.3 - this would have normally been suppressed by the PDF in 
conventional calculations. See Ref. '12D for further discussion about uncertainties 
linked to those approximations. Integrating over Pt, for -^5 = 1.96 TeV with the 
LO MRST2001 gluon PDF'i^ at the scale n = O.^JmI + P^, for |y| < 0.6 and 



Xg = 0.8 GeV, the integrated cross section is 

cr ~ 2.7fib. (42) 

This seems in agreement with the latest measurement by CDfl^ at = 1.96 
TcV in the whole Pt range but for the total cross section only; the extraction of 
the prompt signal was only done for Pt > 1.25 GeV. 

As exposed in Ref. 11201 this calculation is affected by the following uncertainties: 

(1) Choice of the effective gluon mass, A^: for = 0.5, 0.8, 1.0 GeV, a is 2.0, 2.7, 
4.0 nb; 

(2) Choice of the factorisation scale (at which the PDF are evaluated) and 
renormalisation fip scale (at which as is evaluated) : defining /ig = ^ Mq + P^ , 
setting /ii? = /i/j to 0.5 /io, /io, 2/io, cr is 2.7, 2.3, 1.5 

(3) Choice of the cut-off Xmax- its variation introduces NLL corrections in the BFKL 
approach. 

Beside those, we have the usual uncertainties linked to the PDF (especially at 
low x) and a possible iiT-factor or equally higher-order pQCD corrections. 

Using the same parameters and setting y — 0, one can in turn compute the cross 
sections for 'ip{2S), but also for T{nS) at = 1-96 TeV and at = 14 TeV (see 
Table 18). 



Table 18. Direct cross section calculations Ref. I12U1 



da/dy\y^f^ 


J/4> (pb) 


V'(25) (pb) 


T{15) (nb) 


T(2S') (nb) 


T{35) (nb) 


^ = 1.96 TeV 


2.2 


0.6 


40 


12 


9 


^ = 14 TeV 


8.1 


2.5 


310 


100 


80 
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3.4.2. Pt differential cross section 

Unfortunately, one cannot rely on the amplitude written above the compute the Pt 
differential cross section (see Ref. I12L)|| . As a makeshift, they use a simple parame- 
terisation for the partonic cross section based on dimensional counting: 

2^5 ^''(a4"+^t) (43) 



dP^ °^ 5 ^ p2)3 ■ 

Taking again x^ax = 0.3 and normalising the distribution by equating its inte- 
gral over P^ to the previous one, one obtains a reasonable agreement with the data 
from CDF. Again, the comparison is somewhat awkward since their prediction (for 
^/s = 1.96 TeV) is only for direct production (what they call "prompt") whereas 
the data at ^/s = 1.96 TeV are still only for prompt (and for Pt > 1.25 GeV only). 

3.4.3. Other results 

Since some NNLO processes seem to have enhanced contributions, a second class 
of diagrams was considered where the two t-channel gluons "belong to two different 
pomerons" , or two different partons showers. This kind of contributions can be 
related to the single diffractive cross section (see Ref. fl20)) . It was however found 
that this class of diagrams contributes less that the one considered above, though 
it may compete with it at large ^/s. 

In the same fashion, they consider associative production Q + cc, for which they 



expect a{tpcc) to be close to 2 nh, which gives for the hadronic cross sectio''^^ 



da{il:cc) 



= xig{xi)x2g{x2)<yiiljcc) w OM^ib. (44) 

y=o 



dy 

In this case, it is about only 1 % of the other contributions. 

3.5. CES: Comover Enhancement Scenario 
3.5.1. General statements 

In their niodePinSI, P. Hoyer and S. Peigne postulate the existence of a pertur- 
bative rescattering of the heavy-quark pair off a comoving colour field, which arises 
from gluon radiations preceding the hea vy-q uark pair production. The model is first 
developed for low Pt < m production * * and generalised in a second wor 
to large Pt 3> to, where quarkonium production is dominated by fragmentation. In 
general {i.e. at low and large Pt)i the model assumes a rich colour environment in 
the fragmentation region of any (coloured) parton, modelled as a comoving colour 
field. In the case of large- Pt production, the comoving field is assumed to be pro- 
duced by the fragmenting parton DGLAP radiation. The strength and precise shape 
of the comoving field are not required to be the same at small and large Pt. 

If, as they suppose, the presence of the comoving field is responsible for an en- 
hancement of the production cross section, as it is absent in photon-hadron collisions 
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(no colour radiation from the photon and no fragmentation in the Pt region where 
the data are taken) there would not be any increase in photo-production. Indeed, 
the NLO CSM cross sectioiff^ fits well the data and no modification is needed. 

Taking benefit of this assumed perturbative character of the scattering and as- 
suming simple properties of the comoving colour field - namely a classical isotropic™ 
colour field -, they are able to carry out the calculation of the rescattering, even 
when two rescatterings of the heavy-quark pair are required, as in the gluon frag- 
mentation cas^^^, to produce a colour-singlet ■^S'l state. In the latter case, which 
is relevant at large transverse momentum, the tp' (as well as directly produced J/V") 
is predicted to be produced unpolarised. 

Another assumption, motivated by the consideration of the relevant time scales, 
is that the heavy quarks propagate nearly on-shell both before and after the rescat- 
tering. In the latter case, the assumption is comparable to the static approximation 
of the CSM. Furthermore, the rescattered quark pair is projected on a colour-singlet 
state which has the same spin state as the considered Q, similarly to the CSM. 

Note that since the strength of the field is unknown, only cross-section ratios 
and polarisations can be predicted in the fra mew ork of this model, not absolute 
normalisations. Let us review the high-Py cas which interests us most. 



3.5.2. Results of the model 

Considering two perturbative scatterings as illustrated in Fig. 1271 and working in 
a first order approximation for quantities like ^ and ^ for on-shell quarks before 
and after the rescattering, they sho"vA^ that the rescattering amplitude to form 
a colour-singlet ^Si quarkonium from a gluon can be cast in the following simple 
form'^: 

Mi^Su S,) = Z^i£M dab.nGi ■ e{\g) 62 ■ e{S,)\ (45) 
2v Bttto-^ 

where \g is the polarisation of the fragmenting gluon, Gi — Tf{£i) x £i involves 
the colour field Tf produced by the fragmentation, £i and £2 are -as depicted in 
Fig. 1271 - the momenta of the two rescattering gluons. As usual, Ti{0) is the radial 
wave function at origin. Apart from the presence of rescatterings, their approach 
follows the same lines as the CSM with a projection on static colour-singlet states 
with the proper spin state. 

Very interesting pieces of information can be extracted from this formula: since 
Eq. H45|l is independent of F ^ , solely transverse gluons emerging from the comoving 
gluon field contribute to ^Si quarkonium production at high Pr- Gauge invariance 
is also preserved: Ai{T'j{£i) — > £^) — since the quarks are on-shell. Besides this, 

™in the CMS of the heavy-quark pair. 

"the details of the calculation can be found in Ref. IT281 
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the most interesting part is the factorised dependence on the fragmenting gluon 
and quarkonium polarisations. They appear in two different scalar products with 
the quantity G involving the comoving field F f . The polarisation of the quarkonium 
depends solely on the properties of the comoving field through G2- 

This outcome of this p ertu rbative calculation is totally at variance with the 
statement of Cho and WisJS^I j a ■^S'l produced by a high Pt (thus real and 
transverse) gluon is to be transversally polarised. More precisely, if the comoving 
field is supposed to be isotropically distributed in the comoving QQ rest frame, then 
they predict an unpolarised production of '^Si quarkonia in the high Pt region. 




To what concerns P-waves, they predict that "rather independently of the form 
of F/" , xi should be "longitudinally polarised if the rescattering scenario dominates 
Xi production at high Pt" ■ Finally, we should emphasise that the approach at high 
Pt applies equally for charmonium and bottomonium as long as gluon fragmentation 
matters. This sole consideration of gluon fragmentation can be seen as a drawback: 
indeed c-quark fragmentation seems to be bigger than gluon one since it is one order 
of as less (see Fig. □ (left)). 

In conclusion, this si mple production model seems quite interesting (and moti- 
vated by phenomenology! |), even though the properties of this comoving field are 
poorly known (to obtain the results mentioned here, it was supposed to be classical 
and isotropically distributed in the QQ rest frame). The existence of such perturba- 
tive reinteraction should be also searched for in other processes involving quarkonia. 
A qualitative agreement between the predictions of the model and future data on 
polarisation and cross-section ratios would be quite intriguing and might indicate 
the importance of comovers in quarkonium production. 

3.6. Non-static and off-shell contributions to the CSM 

As we have seen, none of the models reviewed above provides an entirely conclu- 
sive solution to the heavy-quarkonium production problem, especially without the 
introduction of new mechanisms and associated free parameters. Considering that 
the CSM is still the most natural way to describe heavy-quarkonium production, 
it is legitimate to wonder whether some of the hypotheses of this model were actu- 
ally justified. One of these is in fact common to the CSM, the COM and the CES 
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and can be referred to as the static approximation, in the sense that one considers 
the hard part of the process solely for configurations where the heavy quarks are 
on-shell and at rest in the frame of the quarkonium to be created. 

This approximation is usually justified by the fact that the wave function, pa- 
rameterising the amplitude of probability for the binding, should be peaked at the 
origin and expanding its product with the hard amplitude, one has only to keep the 
first non- vanishing term in the Taylor series. In the latter, the hard part is evaluated 
at this static and on-shell configuration. 

This is legitimate as long as the hard part is "well-behaved" , e.g. does not present 
singularities in the non-static region. In order to test the validity of this supposition, 
one has to reconsider the basis of quarkonium (Q) production in field theory. In the 
following, we review a consistent and systematic scheme to go beyond this static 
approximation proposed by J. P. Lansberg, J. R. Cudell and Y. L. Kalinovsky. In 
the latter scheme, when off-shell configurations are taken into account, there exist 
new contributions at the leading-order for '^Si production by gluon fusion. These 
can compete with the ones considered in the CSM. 

This approach requires the use of 3-point vertices depending on the relative mo- 
mentum of the constituent quarks of the Q which can be paralleled to Bethe-Salpeter 
vertex functions. The arbitrary normalisation of these is fixed by the calculation of 
the Q leptonic width. Gauge invariance, which could be spoilt by non-local contri- 
butions, is preserved by the introduction of 4-point vertices. 

3.6.1. Description of the bound-state nature of the quarkonia 

The transition qq ^ Q can be described by the following 3-point functior^^I^. 

(p,p) = r(p,p)7^, (46) 

with P = pi — p2 the total momentum of the bound state, and p = {pi + P2)/2 
the relative momentum of the bound quarks. This choice amounts to describing the 
vector meson as a massive photon with a non-local coupling and is justified by other 
studies (see e.g. Ref. ll26|l . 

In order to go beyond the static limit, the quarks are not supposed to be on- 
shell. To allow connections with wave functions, r(p, P) is taken as a function of 
the square of the relative cm. 3-momentum p of the quarks, which can be written 
in a Lorentz invariant form as p^ = —p^ + ^^'j^'i and possible cuts in the vertex 
function T{p, P) are neglected. Two opposite scenarios have been considered: 

N 

r(p,P) = ^ and T{p,P) = Ne-T^, (47) 

(1 + A^') 



both with a normalisati on and a size parameter A, which can be obtained from 
relativistic quark model j ' l -'-^o l J^^g' l J^JUl xhe normalisation N can be fixed from the 
leptonic-decay width and the procedure to do so is thoroughly explained ij ^^ l ^^^ l 
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3.6.2. LO production diagrams 

As usual, it is natural to think of gluon-fusion process to produce quarkonia at 
high-energy, especially if the value of Pt is not excessively large. To conserve C 
parity, in the case of ^Si production, a third gluon is required; this leads to the 
consideration of gg Qg processes. 

As explained in Ref. .12K1 the use of the Landau equations to determine the 
discontinuities of the amplitude leads us to the distinction between two f amilies of 
diagrams, the first is the one which gives the usual LO GSM contributionJ^S^HEi 
when the relative momentum of the heavy quarks entering the Q is set to zero, the 
other family is absent in the on-shell limit and was never considered before Ref. 11251 
in the case of Q inclusive production. These are shown in Fig. [2HI 





^^^^ 




%1 



'■Si , 



\ 




Fig. 28. The first family (left) has 6 diagrams and the second family (right) has 4 diagrams 
contributing to the discontinuity of gg '^Sig at LO in QCD. 



3.6.3. Gauge invariance 

Now, the diagrams of Fig. |2H1 (right) are not gauge invariant. Indeed, the vertex 
function T{p, P) takes different values in diagrams where either the on-shell quark 
or the antiquark touches the Q, so that current conservation is shattered compared 
to the case where the Q is replaced by a photon. 




Fig. 29. Illustration of the necessity of a 4-point vertex. 
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Physically, this can be understood by the lack of some diagrams: if one considers 
a local vertex, then the gluon can only couple to the quarks that enter it. For a 
non-local vertex, it is possible for the gluon to connect to the quark or gluon lines 
inside the vertex, as shown in Fig. These contributions are related to a 4-point 
qqQg vertex, r^1](ci, C2, q, P) , whose form is, in general, unknown, although it must 
obey some general constraint^2SE^33 

• it must restore gauge invariance: its addition to the amplitude must lead to 
current conservation at the gluon vertex; 

• it must obey crossing symmetry (or invariance by C conjugation) which can be 
written 

r^t^ci, C2, q, P, m) = -7or^1](-C2, -ci, g, P, -m)^7o; (48) 

• it must not introduce new singularities absent from the propagators or from 
T{p, P), hence it can only have denominators proportional to (ci — P)^ — or 
(c2 + P)2-to2; 

(3) 

• it must vanish in the case of a local vertex F)^ 7^, hence we multiply it by 

r(2ci-p,P)-r(2c2 + p,P). 

These conditions are all fulfilled by the following simple choic^^^ 

r(fj(ci, C2, P, q) = -igsT^, [r(2ci -PP)- r(2c2 + P, P)] 

' + 1 ^ (m 

_(c2 + P)2-m2 (ci-P)2-m2j '^^ ^ 
where the indices of the colour matrix T arc defined in Fig. I3UI and is the strong 
coupling constant. 




C2 q 



Fig. 30. The gauge-invariance restoring vertex, r^**). 



It should be also emphasised that this choice of vertex is not unique. A first 
study of the effect of different choices can be found in Ref. |^ where we can 
see that agreement with data can be reached. The introduction of such new vertex 
can be paralleled with the COM since the quark pair (ci, C2) that makes the meson 
is now in a colour-octet state. Such configurations are required here to restore gauge 
invariance. 

This new F*^**) vertex introduces two new diagrams in the calculation of the 
amplitude which, now, becomes a gauge-invariant quantity. The detail of the calcu- 
lation of the polarised cross section can be found in Refs. 11251 1811 
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3.6.4. Results for J and ip' 

We show in Fig. EH (left) the results obtained in Ref. [HSI for = 1800 GeV, 
I77I < 0.6, TO = 1.87 GeV and A = 1.8 GeV. The curves for cttot, ctt and ctl are 
calculated within this approach with the new cut only (see Ref. 1 1251 for details), the 
LO GSM is recalculated from the expression of Refs. 1771 1^1^ 

In fact, the normalisation of the results using the decay width removes most 
dependence on the choice of parameters: instead of a factor 100 of difference ex- 
pected from ^ j^^=^ ° j there is less than a factor 2 at Pt ~ 4 GeV and a factor 3 
at Pt = 20 GeV. Interestingly, the dependence on A is negligible once values of the 
order of 1.4 GeV are taken. 




6 8 10 12 14 

Pj (GeV) P,- (GeV) 



Fig. 31. (left) Polarised (or and cr^) and total (cftot) cross sections obtai ned w ith a Gaussian 
vertex function, m = 1.87 GeV, A = 1.8 GeV and the MRST gluon distributioiC2E( to be compared 
with LO GSM. (right); Polarised (or and aj^) and total (utot) cross sections for 1/)' obtained 
with a Gaussian v ertex function. a„ode = 1.334 GeV, m = 1.87 GeV, A = 1.8 GeV a nd the CTEQ 
gluon distribution EH to be compared with LO CSM, GSM fragmentatiorEEl and the data 
from CDlE] 

We see that the contribution of the new cut cannot be neglected at large Pt- It is 
interesting to see that it is much flatter in Pt than the LO CSM, and its polarisation 
is mostly longitudinal. This could have been expected as scalar products of ei, with 
momenta in the loop will give an extra contribution, or equivalently an extra 
Pt power in the amplitude, compared to scalar products involving et- 

Although the use of the leptonic width is a good means to fix most of the 
dependence on N in the J /ip case, it is not so for radially excited states, such as 
the . Indeed, in this case, the vertex function must have a node. We expect that 
^2s{PiP) should be well parametrised by 

1-^^ r-^oriV' 1-^i-J- e^, (50) 

V anode) (1 + Jf^)2 V anode] ^ ' 

where anode is the node position. 

However, the position of the node is not very well-known, and it is not straight- 
forward to relate our vertex with off-shell quarks to an on-shcU non-relativistic wave 
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function. Let us suppose here that the node has the same position in the vertex func- 
tion than in the wave function: now, the integrand changes sign at anode and the 
positive contribution to the integral can be compensated by the negative one for a 
precise values of anode] the latter turns out to be close to the estimated value of the 
node in the wave function. In other words, our normalisation procedure, translated 
from that of the CSM, gives a large value for N compared to the J/t/; case, whereas 
the calculation of the ip' production amplitude is not affected much by the presence 
of a node in T. 

Fig. 1^ (right) shows that for anode — 1-334 GeV, one obtains a good fit to 
CDF data at moderate Pt (note that the slopes are quite similar. This is at odds 
with what is commonly assumed since fragmentation processes -with a typical 1 / 
behaviour- can also describe the data) . The -0' is predicted to be mostly longitudi- 
nal. 

In conclusion, the procedure followed by Lansberg et al. enables to go beyond 
the static approximation of the CSM by taking into account configurations where 
the quarks which form the quarkonium can have a nonzero relative momentum and 
can be off-shell. As one can see for radially excited states like ■0', very important 
effects can be omitted otherwise. Another important point raised is that new cut 
contributions can appear in this non-static extension. These have been computed 
and were shown to be at least non- negligible at high Pt- Indeed, there exist am- 
biguities in the way to preserve gauge-invariance and these might leave room for a 
better description of the data. 

4. Conclusions and perspectives 

In this review, we have tried to provide a wide overview about the problem of 
the hadroproduction of J/i^, 0', T(15'), T{2S) and T(3S'). In order to clarify the 
discussion, we have limited ourselves to the study of direct production, which does 
not involve decays of higher excited states (like P-waves) nor of beauty hadron in 
the charmonium case. 

On the experimental side, we benefit nowadays from studies from two hadron 
colliders, the Tevatron and RHIC. We have reviewed in detail the procedure used 
and the results obtained by the CDF collaboration. The situation is ideal for -0' as 
they are able to extract the direct cross section as well as the polarisation measure- 
ment for the sole direct sample; for J/-0 and T(IS'), they have extracted the direct 
cross section and a polarisation measurement, including excited states contributions 
though; for T(2S'), T(35'), solely the total cross section was measured. To what con- 
cerns the Run II, they are able to go to smaller Pt but the direct signal as well 
as the polarisation measurement is still under study. Preliminary data nevertheless 
tend to show a confirmation that the (prompt) J /ij) are in fact produced unpo- 
larised if not longitudinal. This is in striking contradiction with all the NRQCD 
predictions. To what concerns RHIC collider, the PHENIX collaboration measured 
for the first time the J /ip cross section at y^s — 200 GeV where the non-prompt 
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signal is negligible contrary to Xc feed-down. 

On the theoretical side, many models and calculations have been proposed since 
the outbreak of the ■0' anomaly in the mid-nineties. We have reviewed six of them. 
In order to present them on the same footage, we have limited ourselves to their 
basic points and their most important results. 

We have first started with the Soft-Colour Interaction approach which deals with 
Monte Carlo simulation. It introduces the possibility of colour quantum-number 
exchanges between partons produced during the collision. This effectively opens a 
gluon fragmentation channel into '^Si at LO. The sole new parameter brought in is 
in fact kept at the same value as in the description of rapidity gaps but the approach 
is likely to suffer from some of the drawbacks as the Colour-Evaporation Model. 

Secondly, we have reviewed the NRQCD factorisation approach, which is usu- 
ally embedded, in hadroproduction, in the Colour-Octet Mechanism. We have pre- 
sented the initial motivations of the approach, namely to cure IR divergences of 
the Colour- Singlet Model. Since it introduces unknown non-perturbative parame- 
ters called Long-Distance Matrix Elements (LDME) , we have reproduced the values 
obtained by different groups, as well as the values obtained for the Colour-Singlet 
matrix elements. As said before, this approach is in contradiction with data to what 
concerns the polarisation. Some modifications of the scaling rules, which normally 
provide a hierarchy between the LMDE's, were proposed to reduce the discrepancy 
but the confirmation of the new CDF preliminary measurement would certainly 
blacken the setting. Finally, according to a recent analysis centered on fixed-target 
experiments by F. Maltoni et aL--*^^"^ the universality of these LDME's is far from 
being observed. 

Thirdly, we have reviewed the application of the kr factorisation approach to 
hadroproduction. We have emphasised the two main differences with the collinear 
approaches (CSM or COM), namely the use of unintegrated PDF's and of the 
effective BFKL vertices. As we have said, the LO contributions include naturally 
fragmentation channels; this lets us think that the approach is perhaps more suitable 
than the collinear ones. In the case of COM, we have given the values of the LDME's 
obtained. These are smaller than the corresponding ones fit from the collinear cross 
sectio ns a nd thus in better agreement with what is ex pecte d from ep data from 
HERi^ and perhaps also from fixed-target experiment Jl^. 

Fourthly, we have briefly reviewed a calculation proposed by physicists from 
Durham of some specific NNLO pQCD contributions of the CSM. These can be 
viewed as LO BFKL contributions and thus are expected to be enhanced compared 
to other NNLO pQCD ones. Unfortunately, the method used cannot predict the Pt 
slope of the cross section. 

The fifth model reviewed deals with rcinteractions with comovers. Two pictures 
were analysed: one for small Pt events where the comovers are created by Bremfi- 
trahlung of the colliding gluons and another for high-Py events where they are pro- 
duced by DGLAP radiation of the fragmenting parton. The authors of the model 
suspect such re-scatterings to be source of a possible enhancement of the cross sec- 
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tion, and if the picture holds, they predict that in the case of gluon fragmentation, 
the quarkonia are produced unpolarised. 

The last approach reviewed is a consistent scheme to go beyond the static and 
on-shell approximation used, for instance, in the CSM and the CEM. The loss of 
gauge invariance caused by the introduction of such non-local effects was shown to 
be curable by the introduction of 4-point vertices. In the case of ■0', the presence of 
a node in the vertex function, which is typical of radially excited states and whose 
effects are necessarily neglected in the static approximation, was shown to have 
non-negligible effects and is therefore expected to be analysed in other processes. 

Finally, on a more formal basis, for all models, we are lacking a factorisation 
theorem suitable both for charmonium and bottomonium cases. To what concerns 
NRQCD, even though some modification of it are required , it has been shown re- 
cently that the (modified) picture still holds at NNLC^'^^ at least for fragmentation 
processes which are supposed to dominate in the large Pt regime. 

In view of similar discrepancies between theory and data at i?-factories, we 
strongly believe that the forthcoming efforts both on the experimental side and 
theoretical side will be fruitful and will certainly shed light on the interplay between 
the perturbative QCD and the bound-state physics for which quarkonium physics 
is archetypal. 
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